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ABSTRACT 


Displacement Transfer Functions were formulated for deformed shape predictions of structures under 
large deformations. By discretization of the embedded beam (depth-wise cross section of structure along the 
surface strain-sensing line), the distributions of surface bending strains could be represented with piecewise 
linear functions. Such a piecewise approach enabled piecewise integrations of the curvature equations 
(Eulerian, Lagrangian, and Shifted curvature equations) of the deformed embedded-beam to yield 
closed-form Displacement Transfer Functions, which are expressed in terms of embedded beam geometrical 
parameters and surface strains. By inputting the surface strain data, the Displacement Transfer Functions 
can then transform the surface strains into deflections for mapping out the overall structural deformed 
shapes for visual display. A long tapered cantilever tubular beam was chosen for studying the 
large-deformation shape predictions. The input surface strains were analytically generated from 
finite-element analysis. Also, the finite-element generated slopes and deflections were used as reference 
yardsticks to study the theoretical shape prediction accuracies. The Displacement Transfer Functions based 
on the Shifted curvature equation were found to be amazingly accurate beyond expectation. However, the 
Displacement Transfer Functions based on the Eulerian and Lagrangian curvature equations, resulted in 
very complex mathematical expressions, and gave very poor shape predictions at large deformations. 

NOMENCLATURE 

c depth factor (vertical distance from neutral axis to bottom surface of uniform beam), in 

c(s) depth factor (vertical distance from deformed neutral axis to bottom surface of nonuniform 

embedded beam), in 

c(x ) depth factor (vertical distance from undeformed neutral axis to bottom surface of nonuniform 

embedded beam), in 

C- depth factor (vertical distance from neutral axis to z-th strain-sensing station on bottom 
surface of nonuniform embedded beam), in 
c 0 value of c- at fixed end (beam root), x = x 0 = 0 , in 

c n value of c- at free end (beam Math equationsp), x = x n = / , in 

DLL design limit load 
E Y oung ’ s modulus , lb/ in 2 

i = 0,1,2,3, — ,/? , strain-sensor identification number, or an imaginary unit, i = yf-l 

/ length of embedded beam, in 

n index for the last span-wise strain-sensing station (or number of domains) 

R radius of curvature, in 

R averaged radius of curvature for domain, x M < x < x- , in 
RRF Rotated reference frame 

SPAR Structural Performance And Resizing finite element computer program 
STS Space Transport System for the Space Shuttle flights 
5 curved coordinate along deformed embedded beam elastic curve, in 
u axial displacement, in 

x, y Cartesian coordinates (x in beam axial direction, y in lateral direction), in 

x i axial coordinate of z-th strain sensor (or strain-sensing station at x = x f ), in 
y(x) beam deflection at axial location x, in 

y curved deflection (curved distance traced by the same material point from its initial un- 
deformed position to its final deformed position, in 

y. = y(x i ) , beam deflection at axial location, x = x f , in 
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V =dy/dx = tan 9(x) , beam slope, in/in 

V i = V(x-) , value of V at strain-sensing station, x £ in/in 

M =(x- -x M ) = l / /7 , domain length (distance between two adjacent strain-sensing stations {x M , x-} , 
in 

A/ c chord length of arc formed by bent A/, in 

AL. projected domain length on the x-axis of bent and tilted A/, in 

£( 5 ) surface bending strain at curved locations , in/in 

£(x) surface bending strain at axial location*, in/in 

e, = 4*/), surface bending strain at strain-sensing station, Jt -, in/in 

s s axial strain in 5-direction, in/in 

6{s) beam slope in reference to 5-system, rad or deg 

0(x) beam slope in reference to x-system, rad or deg 

0. = 0(x i ) , slope at strain-sensing station, x = x f , rad or deg 

0 n = 6{x n ) , slope at beam tip strain-sensing station, x = x n = / , rad or deg 

£ = x - x M , local axial coordinate for domain, x M < x < x f , in 

INTRODUCTION 

For structure (for example, aircraft wings) deformed shape predictions using discretely distributed 
surface strains (bending strains), the Displacement Transfer Functions (refs. 1-5) are needed to convert the 
surface strains into out-of-plane deflections so that one can plot the overall structural deformed shape for 
visual display. The surface strains are to be obtained from the strain-sensing stations evenly distributed 
along each strain-sensing line on the surfaces of the structure. The depth-wise cross section of the structure 
along the surface strain-sensing line can then be considered as an imaginary embedded beam. In the 
formulations of earlier Displacement Transfer Functions (refs. 1-5), the embedded beam was first evenly 
divided into multiple small domains, with domain junctures matching the strain-sensing stations. Thus, 
within each small domain, the beam depth factor could be described with a linear function, and the surface 
strain variation could be described with either linear or nonlinear function. Such piecewise approach enabled 
the piecewise integrations of the curvature equation [classical (Eulerian), physical (Lagrangian), and Shifted 
curvature equations] for the elastic curve of each deformed embedded beam to yield closed-form slope and 
deflection equations in recursive forms. Those recursive slope and deflection equations were then combined 
into a single deflection equation in dual summation form (called the Displacement Transfer Function), 
which is expressed in terms of embedded beam geometrical parameters and discretely distributed surface 
strains. By inputting the surface strain data into the Displacement Transfer Function, one can then calculate 
deflections along each embedded beam. By using multiple strain-sensing lines, the overall structural 
deformed shape (induced by bending and torsion) can be geometrically mapped out for visual display. 

The Displacement Transfer Functions (in different mathematical forms for different types of structural 
geometry) developed by Ko (refs. 1-5), can be used to calculate deflections from strain data obtained from 
the fiber optic strain-sensing method introduced by Richards (ref. 6) to create a revolutionary structure- 
shape-sensing technology called; “ Method for Real-Time Structure Shape-Sensing ,” U.S. Patent Number 
7,520,176 (Ko-Richards), (ref. 6). This patented technology is quite attractive to the in-flight deformed 
shape monitoring of highly flexible wings of unmanned flight vehicles by the ground-based pilot for 
maintaining safe flights. In addition, the real time wing shape monitored could then be input to the aircraft 
control system for aero-elastic wing shape control. 

The surface strain data for inputs to the Displacement Transfer Functions (ref. 1-5) for shape 
predictions can be obtained from the conventional strain gages, from fiber optic strain sensors, or 
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analytically calculated from finite-element analysis. It is important to mention that without the use of the 
Displacement Transfer Functions, any type of surface strain sensors can only sense the strains, but not the 
out-of-plane deflections nor the cross-sectional rotations of the structure. 

The distributed fiber optic surface strains can also be input to the Displacement Transfer Functions, the 
Stiffness and Load Transfer Functions to calculate structural stiffness (bending and torsion), and operational 
loads (bending moments, shear loads, and torques) for monitoring the flight- vehicle’s operational loads. 
This patented method is called “ Process for Using Surface Strain Measurements to Obtain Operational 
Loads for Complex Structures ,” U.S. Patent No. 7,715,994 (Richards-Ko), (ref. 7). The accuracy of this 
patented method for estimating operational loads on structures was analytically confirmed by using 
finite-element analysis of different aerospace structures (tapered cantilever tubular beam, depth-tapered 
un-swept wing box, depth-tapered swept wing box, and doubly-tapered generic long-span wing) (ref. 8). 

The earlier Displacement Transfer Function (refs. 1-5) were analytically validated for prediction 
accuracies by finite-element analysis of different sample structures such as cantilever tubular beams 
(uniform, tapered, slightly tapered, and step-wisely tapered), two-point supported tapered tubular beams, flat 
panels, and tapered wing boxes (un-swept and swept). 

To eliminate mathematical indeterminacy at the limit of uniform beam depth, additional First-Order and 
Second-Order displacement transfer functions (ref. 9) were formulated for structural deformed shape 
predictions. The prediction accuracies of the two new Displacement Transfer Functions were found to be 
comparable with those of the earlier Displacement Transfer Functions 
(refs. 1-5). 

Also, by using piecewise nonlinear (instead of piecewise linear) strain representations in the 
formulations of the improved Displacement Transfer Functions (ref. 10), the shape prediction accuracies 
could be improved greatly especially for highly tapered slender structures. 

The earlier Displacement Transfer Functions (refs. 1-5) were formulated for the straight beams. 
However, by introducing empirical curvature-effect correction terms, the resulting modified Displacement 
Transfer Functions could be applied to the shape predictions of slender curved structures (including 
full-circle structures), with great success (ref. 1 1). 

All the past Displacement Transfer Functions were formulated by integrating the beam curvature 
equation in reference to the undeformed fixed coordinate system, and were found to be amazingly accurate 
for the shape predictions of various slender structures under small bending deformations. 

For certain new type of unmanned aircraft such as Global Observer, with highly flexible wing structures 
having long wingspan reaching 175 ft, the wing tip deflection could reach up to 32 ft during flights. For the 
shape predictions of such type of highly flexible structures, there is a need to formulate large-deformation 
Displacement Transfer Functions. 

The present technical publication discusses different curvature equations [classical (Eulerian), physical 
(Lagrangian), and Shifted curvature equations] for the embedded beam elastic curves, referred to deformed 
or undeformed coordinate systems. This technical publication describes the formulations of different 
large-deformations Displacement Transfer Functions through beam discretization and piecewise integrations 
of the said different curvature equations. 

A long tapered cantilever tubular beam was used in the shape prediction accuracy analysis of the newly 
formulated Displacement Transfer Functions. The input surface strains in the present technical publication 
were analytically generated from the finite-element analysis of the tapered cantilever tubular beam. For the 
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shape prediction accuracy studies of the newly formulated Displacement Transfer Functions, the slopes and 
deflections calculated from the finite-element analysis were used as reference yardsticks. 


CURVATURE-STRAIN RELATIONSHIP 

Figure 1 shows the deformed state of a nonuniform beam with changing depth factor, c(s). The 
curvature-strain relationship can be established graphically from figure 1 . The beam elastic curve has local 
radius of curvature, R(s ), within a small beam segment subtended by dO . The un-deformed curve length, 
AB , lies on the beam neutral surface, and the deformed curve length, A'B'{=AB[1+ £( 5 )]} [where s(s) is the 
surface strain], lies on the beam bottom surface. From the two similar slender sectors, O'AB and O'A'B ' , 
the following relationship can be established: 


OA , c(s) 

= 1+-^- 

O'A R(s ) 


AB' 

AB 


= 1 + 4 ?) 


(i) 


From equation (1) the following curvature-strain equation is obtained: 

_L = fV) 

R(s) c(s) 


Equation (2) geometrically relates the local curvature, 1/7?(V), of the deformed beam elastic curve to the 
associated surface strain, £(s), and the beam depth factor, c(s). Equation (2) was the basis for the 
formulations of the earlier Displacement Transfer Functions (ref. 1-5), and is also used in the present 
technical publication to formulate the Displacement Transfer Functions for large deformations. 

DIFFERENT CURVATURE EQUATIONS 

This section discusses three types of curvature equations for deformed beam elastic curve: 1) classical 
(Eulerian) curvature equation, 2) physical (Lagrangian) curvature equation, and 3) Shifted curvature 
equation. These three curvature equations will be geometrically related to surface strains and beam depth 
factor for establishing three curvature-strain differential equations for the formulations of the 
large-deformation displacement transfer functions. 

Classical (Eulerian) Curvature Equation 

The derivation of the classical curvature equation for a plane curve can be found in any standard 
calculus textbook, and are presented in the following paragraphs for the purpose of refreshing memory. 


From figure 1, the infinitesimal curved segment, AB = ds & ^j(dx ) 2 +(dy ) 2 , and the local slope, 

0 = tan ~ x (dy! dx ) , can be used to express the curvature, 1 I R = d6 1 ds , in the following forms (see slender 
sector, O’AB in figure 1) (ref. 12): 


1 _d0_d0dx_y f - ( 1 _ d 2 yidx 2 1 

R(s ) ds dx ds dx [ \dx) dsldx 1 +(dyldxf Jl + (dy / dx) 2 

d 2 y/dx 2 

~{ 1+(</ k /<&) 2 ] 3/2 


( 3 ) 
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Equation (3) is the classical mathematical curvature equation for a plane curve, and can be found in any 
standard calculus textbook. 

When equation (3) is applied to describe the curvature of the deformed-beam elastic curve, the 
x-variable in equation (3) will now denote the abscissa of any material point on the elastic curve in its final 
deformed position (Eulerian description), (fig. 1). Thus, based on equation (3), a certain x-coordinate of a 
material point in the deformed configuration no longer represents the same material point in the undeformed 
state. Namely, the deformed material point has no memory of its initial undeformed x-location. As will be 
seen shortly, the slope equation formulated based on equation (3) gave erroneous slope data, and caused 
deflection predictions to be poor at large deformations. For more detailed discussions of the classical 
curvature equation (3), please see references 13-16. 

It must be mentioned that the classical Euler-Bemoulli beam theory was formulated based on the 
linearized classical (Eulerian) curvature equation (3) [that is, the term (dy I dx) 2 neglected], which was 
referred to the deformed x-coordinate. For small deflections, the difference between the deformed and 
undeformed x-coordinates is negligible; however, for large bending deformations the two x-systems can be 
quite different. 


Physical (Lagrangian) Curvature Equation 

As shown in figure 2, the infinitesimal curved segment, ds, of the deformed-beam elastic curve in the 
limit could be considered as an infinitesimal straight-line segment. Then, the local slope, 0(s ) , along the 

s-curve can be expressed as 0(s) = sin -1 {dy I ds) . The local curvature expressed in reference to the ^-system, 
can then be written as equation (4): 

R(s) ds yds) 


d 2 y/ds 2 


Jl—idvTds^ 


( 4 ) 


In general, under bending deformation, the original undeformed length, dx, along the beam x-axis will 
be deformed into arc length, ds, along the s -curve (fig. 2). Then, one can write: 


ds 

=l + £ 

dx 


( 5 ) 


in which s s is the axial strain in ^-direction. For an in-extensional beam, e(s)=0, and equation (5) becomes: 

ds = dx (6) 

Equation (6) shows that for the in-extensional beam, the .s-coordinatc of any material point in the 
deformed configuration is always equal to its initial undeformed x-coordinate. In view of equation (6), the 
curvature equation (4) can be rewritten in reference to the undeformed x-coordinate as: 

1 d 2 y/dx 2 

R( x ) ^l—(dy/ dx) 2 
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Equation (7) is the physical (Lagrangian) curvature equation for the elastic curve of an in-extensional 
beam expressed in reference to the undeformed x-coordinate (Lagrangian description). Equation (7) was 
also derived by Kopmaz (ref 13) and Hodge (ref 14) based on kinematics of a deformed beam (fig. 2). The 
alternative derivation of equation (7) is described in Appendix A. 

It is important to mention that for an in-extensional beam {ds = dx ), (fig. 2), when the 5-system is 
converted into the undeformed x-system, sin 9{s)(=dy/ds) in the 5-system will become tan 9{x)(=dy/dx) 
in the x-system [that is, sin #( 5 ) = tan#(x) ]. Thus, the slope, 9(x ) , in the x-system will be slightly smaller 
than the slope, 9{s ) in the 5-system. Therefore, the line^'E drawn from point A' (fig 2) to define the 
slope, 9(x ) , in the x-system will be slightly tilted away from the true tangent to the 5-curve at point A' 
(fig. 2). However, this slope-tilt is miniscule because the axial displacements induced by the lateral bending 
are quite small. 


Shifted Curvature Equation 

The formulations of the Shifted Displacement Transfer Functions are based on the undeformed 
x-coordinate. Namely, the deformed material points are shifted back in the x-direction and plotted at their 
respective initial undeformed x-coordinates (fig. 3). Thus, the Shifted curvature equation to be formulated 
will be slightly different from the physical curvature equation (7) for which the deformed configuration of 
the beam is described with the un-deformed fixed coordinate system (that is, no horizontal shifting of 
deformed material points). 

From a small right triangle A A'B'C (fig. 2), the following relationship holds: 

ds 2 = ( dx - du ) 2 + dy 2 (8) 

in which u (= AA" ), (fig 2) is the x-displacement (axial direction) of a material point A. In view of equation 
(8), the physical curvature, equation (4), in the 5-system can be rewritten as: 

1 d 2 ylds 2 d 2 y ds d 2 y ds 

R ( s ) y/l- (dy/ds) 2 ds 2 yfds 2 ~dy dg2 ^{dx-duf 

d 2 y ds d 2 y ds 

= — — = — — (9) 

ds 2 dx —du ds 2 d(x — u) 

Because in the beam bending deformations, the axial displacement, u, of a material point is quite small, 
it is reasonable to set u — » 0 , which is equivalent to horizontal shifting of the deformed material points like 
{A' ,B'} (fig. 3) toward the right to the positions {A" ,B" } with x-coordinates respectively matching the 
initial x-coordinates of the un-deformed positions {A , B} . 


Under the condition, u — > 0 , equation (9) becomes: 

1 _ d 2 y ds 

R(s) ds 2 dx 


( 10 ) 
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For an in-extensional beam, the curved coordinate, s, of any deformed material point is always equal to 
its initial undeformed coordinate x (fig. 3); namely, s = x and ds = dx\ eq. (6)]. Then equation (10) can be 
written in the un-deformed x-system as: 


1 _d 2 y 


( 11 ) 


Equation (1 1) is called the Shifted curvature equation for the in-extensional beam (ds = dx ), referred to the 
undeformed x-coordinate. Equation (11) can also be established graphically in view of figure 3 (see 
Appendix B). 

The mathematical form of equation (1 1) is equivalent to the linearized form of the Lagrangian curvature 
equation (7) [the term (dy/dx) 2 dropped] because both equations (7) and (11) are referred to the 
undeformed x-coordinate (figs. 2, 3). It is important to mention that, for an in-extensional beam, setting 
u —> 0 (that is, shifting deformed material points and plot at their respective undeformed x-locations) causes 

(dy/dx) 2 — >0 according to equation (A8) in Appendix A and, therefore, implying linearization of the 
Lagrangian curvature equation (7). 

However, equation (11) cannot be equivalent to the linearized Eulerian curvature equation (3) [the term 
(dy/dx) 2 removed] even with similar mathematical forms because, equation (11) is referred to the 
undeformed x-coordinate, while equation (3) is referred to the deformed x-coordinate (fig. 1). For small 
bending deformations, the difference between the deformed and un-deformed x-coordinates could be 
negligible; however, for large bending deformations, the two x-systems could be quite different. 

As will be seen shortly, the Displacement Transfer Functions formulated based on the Shifted curvature 
equation (11) turned out to be amazingly accurate beyond expectation for large-deformation shape 
predictions (for example, cantilever beam-tip slope reaching even as large as 66° under bending). 

The unexpected discovery in the present technical publication is that the linearized curvature equation 
used in the formulations of earlier Displacement Transfer Functions (refs. 1-5) was not referred to the 
deformed x-coordinate, but to the undeformed x-coordinate. Therefore, the earlier linearized curvature 
equation was actually the Shifted curvature equation (11), and the Displacement Transfer Functions 
formulated based on equation (11) were found to be amazingly accurate beyond expectation in the structural 
deformed shape predictions. 


CURVATURE-STRAIN DIFFERENTIAL EQUATIONS 

In view of equation (2), the Eulerian, Lagrangian, and the Shifted curvature equations (3), (7), and (11) 
can be written respectively in the forms of curvature-strain differential equations for the deformed beam as: 

1 . Eulerain curvature-strain differential equation: 

d 2 y I dx 2 s(x) 

[i +(dy/<k) 2 ] 3 ' 2 c(x) (12) 
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2. Lagrangian curvature-strain differential equation: 


d 2 y/dx 2 a(jc) 
^(dy/dcf c (*) 

3. Shifted curvature-strain differential equation: 

_ **) 
c(x) 


( 13 ) 


(14) 


The formulations of the large-deformation displacement transfer functions stemmed from the 
integrations of the curvature-strain differential equations (12)— (14). Note that, equations (12)— (14) are 
purely geometrical relationships, containing no material properties. Likewise, the Displacement Transfer 
Functions formulated by integrating equations ( 1 2)— ( 1 4) will also contain no material properties. In fact, the 
material properties will affect the outputs of the surface strains, s(x) (refs. 1-5). 


NONUNIFORM AND UNIFORM BEAMS 

In the present technical publication, the term “beam” implies the imaginary embedded beam (the 
depth-wise cross section of the structure along the surface strain-sensing line), and not the classical isolated 
Euler-Bernoulli beams. The nonuniform beam in the present technical publication is defined as the 
embedded beam with a varying depth factor, c(x ) , a vertical distance from the neutral axis to the bottom 
surface of the nonuniform beam at axial location, x. For a nonuniform tubular beam, c(x) will be the 
outward radius. The uniform beam in the present technical publication is defined as the embedded beam 
with constant depth factor [that is, c(x) = c ]. 

SHIFTED DISPLACEMENT TRANSFER FUNCTIONS 

In the formulations of earlier Displacement Transfer Functions (refs. 1-5) for structural deformed shape 
predictions, the curvature-strain differential equation used turned out to be equivalent to the 
Shifted-curvature-strain differential equation (14). Therefore, there is no need to reformulate the 
large-deformation Displacement Transfer Functions based on the Shifted curvature-strain differential 
equation (14). The Displacement Transfer Functions formulated earlier can be called Shifted Displacement 
Transfer Functions, which, as will be shown shortly, were found to be amazingly accurate beyond 
expectations for the shape predictions not only for small deformations, but also for very large deformations. 

The Shifted Displacement Transfer Functions formulated earlier for equal domain length ( A/=constant) 
(refs. 1-5) are reproduced in the following sections. Keep in mind that the slopes, 0 t [= 0(x t ) = (dy / dx) t \ 
(/ = 1,2,3, — ,/?) , at the strain-sensing stations, x = x f , appearing in the following Shifted Displacement 
Transfer Functions are referred to the initial undeformed x-coordinates. 

Based on Piecewise Linear Strain Representations 

The following four sets of Shifted slope and deflection equations [eqs. (1 5)— (1 8)] were formulated by 
using piecewise-linear functions to describe the distributions of both surface strains and depth factors, 
{<x),c(x)}. 
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1. Nonuniform Beams (refs. 1, 3) 


Slope equation: 


tan 9 • = A / 


'i - 1 


-i-\ 


£i-l C i 


(Cm-c,) 2 


-log— i- 


+ tan 9- 


C, 


i - 1 


(z' = l,2,3,....,zz) 


M 


(15a) 


Deflection equation: 
a. In recursive form: 


y ,=( a /) 2 


2 ( c /-i - c /) (c M -c ; ) 3 


C; log— + (C,_i -C ; ) 


"Z-l 


+ + A / tan 9^ 


i - 1 


(15b) 


(/ = 1 , 2 , 3 ,. ...,«) 


b. In dual summation form: 


= 


= ( 4 ') ! Z 


7-1 


Sj-lCj-SjCj-l 


^ 2(c, 1 -c,) (c,,-c,) 


c z lo g e 


C 7-l 


+ ( C /-1 ~ c /) 


/-I 

\2X^ 


+ (A/) Z > (z-/) 


i=i 


Contributions from deflection terms 


! -£■ £■ i C, - £,C, i C 7 - 

7~‘ 7 _|_ J J j^g J_ 




(c-t-c-r 


■V-l 


> + jv’q + (z)A/ tan 9q 

=0 for cantilever beams 


(15c) 


Contribution from slope terms 


(z = l,2,3,....,«) 


Equation (15c) is called the Nonuniform Displacement Transfer Function. 

2. Slightly Nonuniform Beams [log e terms in equations (15a)-(15c) expanded] (refs. 1, 3) 

Slope equation: 


tan 0. 

Deflection equation: 
a. In recursive form: 


A/ 


2c ; _i 


r \ 
2 — — 
C i~ 1 J 


*/- i + ^ t 


-i-tan ; (z = 1,2,3,. ...,«) 


vz 


(A /) 2 


6c 


7-1 


Si - 1 + ^7 


u 7-i y 


+ + A/ tan ; (z = 1,2,3,. ...,«) 


(16a) 


(16b) 
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b. In dual summation form: 



fi 

C ■ , 

]}- “ ) 
C ■ , 


} +<a ' ) 2 f| 

^ 7 =1 

f O' - 7) 

C ■ , 

Ml 



7-1 

V y- 1 y 



7-1 

_V 7-1 ^ 



Contribution from deflection terms Contributions from slope terms 

+ y 0 + Al tan 0 O 

=0 for cantilever beams 


(16c) 


(/ = 1,2,3,...., «) 


Equation (16c) is called the Slightly-Nonuniform Displacement Transfer Function. 
3. First Order Expansion of 1 / c(x) (ref. 8) 

Slope equation: 


tan ft 


A / 


' 6 c M 


c ; 


c i - 1 y 


s i - 1 + 


5-2-^- 

0-1 


+ tan ; (/ =1,2,3,....,/*) 


Deflection equation: 
a. In recursive form: 


JE = 


(A/) 2 


12c, 


Z-l 


A \ f 

5- — 


O-i y 


O-i + 


3--^- 
O-i y 






+ j,_j + A/ tan ; (z = 1,2, 3,. ...,«) 


(17a) 


(17b) 


b. In dual summation form: 

(A /) 2 v-i (A /) 2 


Tz =- 


12 M C J - 1 


A A 

5- — 

C-i; 


+ 


A 

3 - — 

C-i 7 


Contributions from deflection terms 


+ ■ 


1 

1 

<N 

r 4 _ s- 1 

f /-i + 

f 5 _ 2 S ' 

^ c 

y=i c i - 1 

i 

l C--iy 


7 


=0 for cantitever beams 


Contributions from slope terms 


+ v 0 + A/tan6> 0 (17c) 

(z = 1,2,3,.. ..,«) 


Equation (17c) is called the First-Order Displacement Transfer Function. 
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4. Second Order Expansion of 1 / c(x) (ref. 8) 


Slope equation: 


tan 0- 


M 


1 12 c 


i - 1 


9-4-^- + 

C i - 1 


f \ 
C , 


v c m y 


^•-i + 



f \ 

2 


10-^- + 3 

Ci 


£ i > 

C;- 1 

V C i-1 j 

i 



+ tan# M (18a) 

(* = 1,2,3,. ...,«) 


Deflection equation: 
a. In recursive form: 




(A /) 2 

60c, •_! 


27 -9-^- + 2 


7 a 2 
c ; 


"/-l y 


e i-\ + 



f \ 

2 

> 

11-^- + 3 

Ci 


£ i 

C i-1 

V Ci-1 j 




+ + A/ tan 0 M ( 1 8b) 

(z = 1,2,3,.. ..,«) 


b. In dual summation form: 


jy = 


(A/) 2 1 

60 


27-9— ^- + 2 


c ;-i 


7 A 2 

_2l 

V C 7-V 


f y-i + 


18-11 — — + 3 


f \ 

2 


C ■ 



J 


£ j > 

If /- 1 J 


J 


, (A /) 2 y 1 (z ~ j) 
12 £ C-i 


Contributions from deflection terms 
C 


9-4 — — + 

C H 


\ 2 
c j 




£ j - 1 + 


13 — 10— — h 3 
< 7-1 


vC-i; 


Contributions from slope terms 


+ )?q + (z)A/ tan Oq 


=0 for cantilever beams 


(18c) 


(z = 1,2, 3,. ...,«) 


Equation (18c) is called the Second-Order Displacement Transfer Function. 

Based on Piecewise Nonlinear Strain Representations 

The following three sets of Shifted slope and deflection equations [eqs. (19)— (21)] were formulated by 
using piecewise-nonlinear function to describe the actual distribution of the surface strain, s(x) , but using 
piecewise-linear function to describe the actual distribution of the depth factor, c(x) . 
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1 . Improved Case (ref. 9) 
Slope equation: 


\l Q 

tan 0 , = — — [( 2 c i ~c i _ i }(c i e i _ 1 — 2 c i _ 1 s i ) + cp^s^ log e 

c i-i 


' 2 (c-c^f 

A/ 

Mc-c^f 


[( 5 c, 3 c^,) 6 ^_j 2 ( 3 c t + (c, + c M )^ +1 ] + tan 6 t_, ( 19 a) 


Deflection Equations: 
a. In recursive form: 


T; = 


(A/) 2 


2 (c,-Cm)‘ 


-[(2c,. - c M )(c,c M - 2c ; _ l c i ) + c,.c M ^. +1 1 c,. log,, -7_ - (c,. - c M ) 


(A/) 2 


12(c,- — c M )‘ 


-[(8c,. - 5c, - 2(5c,. - 2c, )e i + (2c,. + c M )£, +1 ] + + A/ tan<9 M 


( 19 b) 


(/ = 1,2, 3, ••••>«) 


b. In dual summation form: 


y,=(A0 2 X 

j = 1 


^ [(2c - Cy_! )(c£ - 2 Cj^Sj ) + c jCj _p j+x Cj log, - (c, - Cj _! ) 

2(c / -c, 1 ) 


C-i 


12 C c 7“ c t-i) 


y t 8 C - 5 C'-i ) £ j - 1 - 2 ( 5c y - 2c /-i 7 + ( 2c j + C-i )^' + i 1 


Contribution from deflection terms 


i-1 


l -—r be - e-i Xee-i - 2c 7-ie ) + c j c j-i £ m ] lo §i- 


+ (A/) 2 ^(i-y> 

7=1 


+ y 0 + (z)A/ tan (9 0 


=0 for cantilever beams 


~7-l 


2 (e-e-i) 

- ~ 7 } r — 7 be - 3 e-i )e-i - 2 ( 3 e - e-i )e + (e + e-i >e + i 1 

^\ c j c j- 1 1 

Contributions from slope terms 


( 19 c) 


(/ = 1,2,3,....,/?) 


Equation (19c) is called the Improved Displacement Transfer Function. 


12 


2. Log-expanded Case [log e terms in equation (19a, b, c) expanded] (ref. 9) 


Slope equation: 


tan#. 


/V 


12c 


i-l 


{[ 5c h + 4(c i -c i _ I ) 2 ]f i _ 1 - Sc i _ 1 (c i -2c i _ 1 )s i +c i _ 1 (2c i - 3c H )er w } + tan#,., 

O' = 1,2,3,...., «) 


Deflection equation: 
a. In recursive form: 


(A If \[lcU + (c, -c m )(8 C< c m -3c, 2 )-3( C< 3 | 

24c U [+ 2c ( ._, [(3c 2 + {c i -c m )(3c,. - 8c w )]c ; - c M [cf + 2(c, - c M ) 2 ]c, +1 J 
+ J’,-i + A/tan# M 

O' = 1,2,3,.. ..,«) 


b. In dual summation form: 


1 ■ 
24 

+ 2 c >-i | 

+(c, - Cj-I )(s<v<v-, - 3cj) - 3(, 
[(3<4+(c J -c J . 1 X3c J -fc >1 )] 

3 3 ^ 

'/ - 6 V-r). 

K- C J-'[ 

I s /-. 

c;+2(c / -c / .,) ! ] 

> 

12 ^ c 3 
/=! C /-l 

Contributions from deflection terms 

[5Cy_, +4(C, ~Cj _ ,) j ~ 8Cy_l( c y — 2C/-l 

+Cy_l(2Cy — 3Cy_,)f /+ , 

)s, 

y f+ y 0 +0‘)A/tan# o 

v v ' 

J =0 for cantilever beams 


Contributions from slope terms 

0 = 1,2,3, — ,«) 


Equation (20c) is called the Log-Expanded Displacement Transfer Function. 
3. Depth-expanded Case [expansion of l/c(x) ] (ref 9) 

Slope equation: 


tan 0: = 


A/ 

24c 


i-l 





r \ 

2(5c, 

1 + 

1 

+ 

1 

oo 

-3--1 



V C i-l J 

f 

A 

2 


C i 

1 

(c M - 12c, 

+ 

+ 

1 

Um 

) 



(3c,-, + 10c, - s M ) 


+ tan 0; 


i-l 


0 = 1,2,3,.. ..,«) 


(20a) 


(20b) 


(20c) 


(21a) 
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Deflection equation: 
a. In recursive form: 


yt 


(A/) 2 


24c,-, 




V C i-1 


1 

+ — 
5 


b. In dual summation form: 


C i-l ) 


12£; 3e t+l ) 


(4^,-1 + 7f; -£,- +1 ) 
+ v,_i + A/ tan 


(21b) 


(/ = 1 , 2 , 3 ,. ...,«) 


A; 


(A/) 2 y-i 1 
24 


(7^+6^.-^) 


-^-1 


(4^+7^-^,) 


1 

+ — 
5 




1 

V c i-i 7 


V 7- 1 

+ 12 *,-- 3*, - +1 ) 


, (A/) 2 y(i-j) 
24 £ Vi 


Contributions from deflection terms 


^ A 

^-1 


2(5s h +8s j -£ j+1 ) 


V C J - 1 7 


(3f H +10f 7 -f /+1 ) 


V.C-1 7 


(*H -12^+3^,) 


Contributions from slope terms 


+ j7q + (z)A/ tan 0 q 

=0 for cantilever beams 


(21c) 


0’ = l,2,3,....,w) 


Equation (21c) is called the Depth-Expanded Displacement Transfer Function. 

Characteristics of Displacement Transfer Functions 

In each of the Displacement Transfer Functions [eqs. (15c)-(21c)], the deflection, at the 

strain-sensing station, x., is expressed in terms of the inboard beam depth factors, (v„ C 2 , ,c t ), and 

the associated inboard surface strains, (<c 0 ,c;,,c; 2 , , ) . including the values of {c t , £ t } at the 

strain-sensing station X i where deflection, y f , is calculated. It is important to mention that equations 
(15c)-(21c) are purely geometrical relationships, containing no material properties. In fact, the values of the 
surface strains, can be affected by the material properties and internal structural configurations. 

Thus, in using equations (15c)-(21c) for shape predictions of complex structures such as long-span 
aircraft wings, there is no need to know the material properties, nor the complex geometries of the internal 
structures. As will be seen shortly, the Shifted Displacement Transfer Functions [eqs. (15c)-(21c)] turned 
out to be amazingly accurate beyond expectations for the shape predictions of a long cantilever tapered 
tubular beam under large bending deformations with beam tip slope even reaching as large as 0 n ~ 66 deg. 
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EULERIAN FORMULATION OF DISPLACEMENT TRANSFER 

FUNCTIONS 


The next attempt is to formulate the large-deformations displacement transfer functions based on the 
Eulerian curvature-strain differential equation (12). Because equation (12) is referred to the deformed 
coordinate system (Eulerian description), the x-coordinate used in this section must be considered as the 
deformed x-coordinate (no change in y-coordinate). It is natural to expect that by retaining the term 

(dy/dx ) 2 , the shape prediction accuracy could be improved; however, as will be seen shortly, the results 
turned out to be entirely negative at large deformations. 

Projected Domain Length 

Under large deformations, strain-sensing stations, x-(z = 1,2,3,....,/?) , on the un-deformed configuration 
will move along the x-axis to their deformed positions. Therefore, in the present Eulerian formulation, one 
must keep in mind that X t now denotes the deformed locations of the strain-sensing stations on the x-axis. 

Similar to the earlier formulations of the Displacement Transfer Functions (refs. 1-5), the undeformed 
beam was first divided evenly into n number of small domains (domain length, A / = //«), so that under 
deformation, the surface strains and depth factors, (e(x), c(x)}, within each deformed small domain, 
x-_j <x<x n could be represented with linear functions (piecewise linear assumptions) of the following 
forms: 


< x ) = £ i - i -( £ i - 1 ~ s i ) - A ^ i ~ 1 ; (*,1 < * < *,) ( 22 ) 

c ( x ) = C i-1 - (Pi- 1 -Ci) - A ^ il ; (*,_! <*<*,) (23) 

In equations (22) and (23), M } and are respectively the values of {^(x),c(x)} at the 

deformed strain-sensing stations, {x t _ x , x-} ; and A L z (= x t - x i _ l ) is the projected domain length along the x- 
axis of the bent and tilted domain length, A / (unchanged for in-extensional beam), (fig. 2), and can be 
calculated from equation (24) (see Appendix C): 

(24) 

(i = 1,2,3,....,/?) 

in which 0 j _ l [= tan ~\dy / dx)^ is the slope at the deformed strain-sensing station, X t l (fig. 4) and can be 
calculated from the yet-to-be formulated slope equations (29) or (34) in the following section. 

The piecewise linear assumptions [eqs. (22) and (23)] enabled the piecewise integrations of the Eulerian 
curvature-strain differential equation (12) to yield the slope and deflection equations for each small domain. 


AL- = 2 


■"/-l 


+ c 7 . 


\ £ i - 1 + £ i J 


sin 


^ A / £j_ x + £j ^ 

2 c M + c- 


cos 


6>._i + 


A/ g M + €j 
2 c M + c x . 
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Slope Equation 

Let V = dy/dx = tan 9(x ) , then the Eulerian curvature-strain differential equation (12) can be written as: 


dV / dx _ f(x) 

[, + ^]“ = 40 <25) 

In view of piecewise linear assumptions [eqs. (22) and (23)], the piecewise-integration of equation (25) 
within each deformed small deformed domain, x t _ x < X < Jt /? can be written as: 


c y dV 


r 

J */i 


S l - X - i £ l - X ~ £ i ) 


X-X, 


i - 1 


AL, 


Ci-l-iPi-X-Ci) 


x - x t _ 

az: 


-dx 


(x H £xSx f ) 


(26) 


in which the lower limit of integral, V t _ x , is the value of V at the deformed strain-sensing station, x f _ 1? 
namely, V t _ x = V (x M ) = tan 6{x i _ x ) = tan 6 lX . 


Integration of equation (26) to obtain the slope equation in closed form is relatively easy; however, 
subsequent integration of the slope equation to obtain the deflection equation in neat closed-form turned out 
to be a nightmare, and an alternative simplified approach had to be introduced as discussed in the 
subsequent sections. 

After carrying out integration of equation (26), one obtains equation (27) (ref 17, and Appendix D): 


V 

V l + F 2 


l 


V LL 
l + K 


- + - 


i-\ 


'i - 1 


"i-\ 


(x-x i _ l ) + M i 


■ S:C, 


&i-l 


s i-\ c i 


( C i ~ C i -\ ) 


log. 


' 7-1 


+ 1 


V C /-1 


A L t J 


(27) 


By squaring both sides, rearranging, then taking the square root, and writing V (x) = tan 0( x ) , one obtains 
the slope equation in the following form (see Appendix D): 


tan 0(x) 


tan 6, 


■yjl + tan 2 ~d~_ x 


■-M — (x - x M ) + AL ; 


( c i - c ;-i ) 


log £ 


( ^ 

5l-1 


X - X j 

c i - 1 ) 


;-l 


+ 1 


1- 


tan 6, 


M (x - x,_i ) + AL, £, ~ lCl log, 


1 + tail 2 0 j{ C i C i- 


( c i ~ c i-l ) 


- S --1 

v c ;-i y 


x - x, 


i - 1 


AL,. 


+ 1 


(28) 


0,--i ^ x < x ; ) 


in which, when / = 1 , for a cantilever beam we have tan 6 X _ X = tan 9 0 = 0 at the fixed end. 
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At, the deformed strain-sensing station, we have x- - x M = A L • , tanft^) = tan ft, , and equation (28) 
becomes: 


tan ft 


M 


tan ft, = 


^1 + tan 2 ft M 


+ AL; 


£ i £ i- 1 _|_ A-l C / A C z-l j 0 g 
C z _ C z -1 ( C z — C z- 1 ) C z -1 


tan ft. 


u 


i - 1 


+ AL ; 


+ tan ft 


i-i 


£ i ft'-l + g i-l C i g j c i-l | 0 g C / 




C, 


/— 1 


(29) 


(/ = 1,2, 3, 


Equation (29) is the Eulerian slope equation for large deformations of nonuniform beams, but is not 
applicable to uniform beams (c, = r, , = c) , because the logarithmic terms, log ( , (c, / c; , ) , and the 
denominators containing factor, (c i - c M ) , will cause mathematical indeterminacy (that is, 0/0). The way to 
avoid such mathematical indeterminacy is to expand the logarithmic function in series form as described in 
the subsequent sections. 

At, the deformed strain-sensing station, x i h we have x ( _, - x ( _, = 0 , tan ft(x,_, ) = tan ft_, , then 
equation (28) becomes an identity equation [eq. (30)] as: 


tan ft 


M 


tan ft_j = ■ 


^/l + tan 2 ft_j 


= tan ft, 


i-i 


(30) 


tan ft 


i-i 


^/l + tan 2 ft_j 


tanft=0 ; (/ = 1,2,3,. ...,«) 


which confirms the mathematical accuracy of equation (28). 


Deflection Equation 

The deflection equation can be obtained by integrating the slope equation (28) as: 


><*)=( tan 0{xyk+y t _ x 


tan ft . £ { — . . . r £. .a - ^ c . , 

I y—+ ( x ~x i _ l )+AL i 11 ‘‘ l log ( 

y 1 + tan ft , c,-C|-i (ft-fi,) 


r c, ,^1 x - x- , 

I ' - 1 1 11 + 1 

U-, j AT, 


-ftx: 


'*ii 


tanft 


i-i 


7l+tan 2 ft, ft-c,, 


(x-x^+Aft ^ ^ log, 




r c. ^x-x,. 


Ur‘J 


+1 


+ T,! 


AT. 


(*, < x < x fl ) 


(31) 
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in which y t _ x = y(x t _ x ) is the deflection at the deformed strain-sensing station, x f-1 . When i = 1 , we have 
y x _ x = y 0 = 0 at the fixed end of a cantilever beam. 

The powerful Wolfram Mathematica Online Integrator (http://integrals.wolfram.com/index.jsp) 
indicated that no integrated formula existed for the mathematical form of equation (31), which contains 
logarithmic functions. 


Log-Expanded Slope Equation 

As indicated earlier that in the slope equation (28), the logarithmic terms and the denominators 
containing factor, (c. - c i _ ] ) , will cause mathematical indeterminacy (that is, 0/0) for the uniform beam case 
(c- = c i _ l = c ) . The way to avoid such a mathematical breakdown problem is to expand the logarithmic 
function in series form as described below. 


Note that in equation (28), the two conditions, 0< 


x - X- 


i-\ 


A L 


< 1 and - 1 < 


f \ 

-5—1 

V C M J 


< 0 , cause the 


required condition, - 1 < 


V c ;-i j 


x - x. 


i - 1 


AL, 


<1, to hold for series expansion of the logarithmic term in 


equation (28) in terms of (c i - c t _ x ) (ref. 16 and Appendix D) as: 

i°g e 


f \ 




-1 


V C M 


x - X 


i-\ 


A L 


+ 1 


- c. A x - x 


i - 1 


-i - 1 


AZ, 


/ \ 
c / -c M 

2 

C \ 

v C <-1 y 


l ^ J 


+ . 


(32) 


Substitution of equations (32) into the slope equation (31), causing the factor, (c- - c -_ } ) , to be canceled out 
to yield (see Appendix D): 


tan 0(x) = 


4 


tan 0, , | g M (j . x ^ (x - x M ) 2 


1 + tan 2 0,_ j c i- i 


-i-i 


2AL, 


tan 0 M | g M (x x ^ Sj-iCt-SjC^ (x-x M ) 2 


1- 


-^1 + tan 2 0 M c i-i 


-i - 1 


2AL, 


(33) 


(*/-l < X < X f ) 


At the deformed strain-sensing station, x = x- , we have x- - x M = AL Z , and equation (33) takes on the form 
(see Appendix D): 


tan + A L, 


tan 0, = ■ 


tan 2 0 M 2 c m 




2 *- 

Si - 1 + Si 

_V C i-1 J 



1 -- 


tan 0,_i 




: + - 


AL, 


1 + tan 2 0 f _j 2c, -i 


2 -- 




g<-i + £ i 


(34) 


(/ = 1 , 2 , 3 , ....,«) 
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Equation (34) is the log-expanded Eulerian slope equation for large deformations of nonuniform beams, 
including uniform beams as special cases because the logarithmic functions and the factor, (c t _ x - c t ) in the 
denominator of equation (28) have been eliminated. 


Log-Expanded Deflection Equation 

The log-expanded deflection equation can be obtained by the integration of the log-expanded slope 
equation (33). 

Let A, B, C, and ^ be defined respectively in equation (35) as: 


, _ tan#,-, 

sfl + tan 2 8, 


B = ^ 


i - 1 


c i-i 


c = — 


2AL, 




(35) 


then equation (33) can be written in a compact form as: 

A + + C<f 


tan 6{x) = ■ 


V i-04 + ^ + af) 2 

The deflection, y(x), can be obtained by integrating equation (36) as: 


f x s\s A-\- + C£ 2 

X-0 = f tan 0(x)dx = f , d% +y i _ 1 

Jx '' Jo Vl -(A + B^+Cff 


(36) 


(37) 


The powerful Wolfram Mathematica Online Integrator (http://integrals.wolfram.com/index.jsp) can 
now carry out the integration of equation (37) to yield the log-expanded Eulerian deflection equation. 
However, the resulting integrated formula is nearly ten pages long, containing various mathematical 
functions including inverse sine functions, elliptic integrals of the first, the second, and the third kinds 
(ref. 17). Because a very accurate alternative simplified deflection equation could be introduced (see the 
following section), no attempt was made to program and use the complicated ten-page-long integrated 
formula of the log-expanded Eulerian deflection equation (37). 


Simplified Deflection Equation 

To bypass direct integration of equation (37), one can introduce a simplified deflection equation by an 
alternative approach with high accuracy. Because the slope, tan#(x), changes very slowly within each 
deformed small domain, x t _ x < x < x i? it is reasonable to represent tan#(x) with a linear function as: 


tan 9{x) « tan O iX + (tan 0 t - tan O iX ) - * l - ; (x fl <x< x x ) (3 8) 

AL t 

in which tan6t can be obtained from the exact slope equation (29) (with logarithmic terms) or from the 
log-expanded slope equation (34). 
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Equation (38) can now be easily integrated to yield the following closed-form alternative deflection 
equation: 


Xx)=f X tan 0(x)dx+y i _ l 

J *i-1 

=tang i _i(x - x^) +(tan6>- -tanfl^) — 2 ^ l) +y i _ l ; (x M ^x^x f ) (39) 

in which y t _ x = ^(x^) is the deflection at the deformed strain-sensing station, x M . 

At the deformed strain-sensing station, x- , we have x i -x i _ l = AL / , and equation (39) gives the 
deflection, j;(jq ) = y i9 in the following simple form: 

= ^-(tan 0 M + tan 0 t )+ y,_, ; (/ = 1,2,3,....,/?) (40) 

Equation (40) is the simplified Eulerian deflection equation for the large-deformation shape predictions, and 
turned out to be very accurate. The high accuracy of equation (40) was verified in the following way. 

Using the slope data calculated from the Shifted slope equation (15a) (very accurate) as inputs to 
equation (40), withAL. replaced with A/, the calculated beam-tip deflection of a tapered cantilever beam 
(table 1) was found to be only 0.0277 percent off from the deflections calculated from the Shifted 
Displacement Transfer Function [eq. (15c)] (very accurate). The high accuracy of the simplified deflection 
equation (40) thus eliminated the need to program and use the ten-page long integrated formula of equation 
(37). 


Table 1. Geometries of the long tapered cantilever tubular beam. 


/, in. 

(Length) 

t , in. 
(Wall 
thickness) 

c 0 , in. 

(Root depth 
factor) 

c m in. 
(Tip depth 
factor) 

Cyi / Co 

(Depth 

ratio) 

a{= tan 1 [(c 0 - c n )/l]}, deg 
(Taper angle) 

300 

0.02296 

4 

i 

1/4 

0.5729 


L AGRAN GIAN FORMULATION OF DISPLACEMENT 
TRANSFER FUNCTIONS 

Lastly, the Lagrangian formulation of the displacement transfer function for large deformations is based 
on the Lagrangian curvature-strain differential equation (13), which is referred to the undeformed 

x-coordinate. Like Eulerian formulation, by retaining the term (dy / dx) 2 in equation (13), one will 
naturally expect more accurate shape predictions; however, as will be seen in this report, the results turned 
out to be very discouraging. 
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Slope Equation 

Writing V = dy l dx = tan 6(x ) (fig. 2), the Lagrangian curvature-strain differential equation (13) can be 
written as: 


dV / dx _ s(x) 

Vl-E 2 c(x) (4 

In light of the piecewise linear assumption of the surface strain and the depth factor, {^(x) , c(x)} [ A L t 
in eqs. {(22), (23)} replaced with A/], integration of equation (41) within the small domain, x i _ l < x < x f , 
can be written as: 


x-x 


i i 


c v dV r x £ ‘ l £ i* A/ 

f -A T =J ; (*,_,<*<*,) 

J ^-i V 1 — Jx 


2 x X - X ( _, 


(42) 


C iX~( C iX- C i) 


M 


On the left hand side of equation (42), the lower limit of integration, V t _ x , is the value of V at the 
strain-sensing station, x l _ l [that is, V i _ l = (dy / dx) x i = tan#(x M ) = tan 9 i _ l ]. 

Carrying out integrations of both sides of equation (42), one obtains the slope equation [eq. (43)] in the 
following form (see ref. 1 and Appendix E): 


. -1 T/ • -1 JT 1 — / x A 7 £i-\ c i ~ C i-\^i 

sin V = sin V t _ x -\ — — — (x - x M ) + A / 


log 


( c /-i ~ c i) ’ ,_1 ' 

from which the slope, tan#(x)[= F] can be expressed explicitly as: 


( c m ~ c -) 
A/c,_i 


(x-x w ) + l 


(43) 


(x,-_i < X < Xj) 


tan 0(x) = sin< 


( c i~ i - c i) 


(x-x M ) 


+ A/ g '" |C ' CM f f log 

(cm-c,-) 2 


(g,--l -c,-) 

A/c ; _i 


(x-x M ) + l 


-i-sin 1 (tan ) 


(44) 


(x f l <X<X,) 

At the strain-sensing station, x t , we have x, - x M = A/ , tan #(x ; ) = tan 0 t , and equation (44) becomes: 


tan Q. = sin< A/ 


£ <-l £ i | g i-l c > c i-l g i [ Q „ 


.( C M _C l) ( C M- C i) C / 


/-I 


+ sin '(tan # M )} 


(45) 


(i = 1,2,3 


As will be seen in this technical publication, the values of the slopes, <9, , calculated from equation (45) did 
not turn out to be as accurate as the corresponding slopes calculated from the Shifted slope equation (15a) 
for large deformations. 
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Deflection Equation 

The deflection equation can be obtained by integrating the slope equation (44) as: 


y(*)= f tan d(x) A 

J */-l 

sin _, (tan6 , i _ 1 )+ £ ‘ l S ‘ 


-r 


sin 


(*-*/-!> 


£ l-l c l~ c l-l £ i 


+ AI 




log 


c, I ~ c i) 
Ala , 


(*-*U> +1 


dx+K-i (46) 


in which y t - 1 is the deflection at the strain-sensing station, X t _ v 

Let A , B , C , D , and £ be defined respectively in equation (47) as: 


; C = A/ CM f f ; C;) 


J = sin 1 (tan 6 , ; _, ) ;B = 


A /c M 


(c f -i - c, ) ' (c M -c,) 

then equation (46) can be written in a compact form as: 

><*)= f^sinp + Clog(D^+ 1)]^+^ ; (0 < £ < A/) 


; £ s x - x M (47) 


(48) 


Using the Wolfram Mathematica Online Integrator (http://integrals.wolfram.com/index.jsp), equation (48) 
can be integrated to yield the following complex form of the Lagrangian deflection equation: 


y(Z)=-L(DZ + iy iC 

In 


iB(Dg + l) 


-iC 


(BD^ + B) 2 


D A 


-i iC 


(- B Y 


r 

\ 


COS 

A -=■ 

- i sin 

A- = 



l D ). 


l D) 


x -r 


D 


Incomplete gamma function 


~(D4 + 1) 


2/C 



f— B 


f 

B V 

" 

COS 


+ i sin 

A 


x r 


l D) 


V 


_ 


D 


+ y t ~ i 


(49) 


Incomplete gamma function 


(0 < £ < A/) 


In equation (49), except for the term y._ x (z is an integer), the symbol, z, appearing in the rest of the 

terms is not an integer, but an imaginary unit, i = a/-T . Equation (49) contains imaginary unit, i = V-T , 

f 00 i z 

and incomplete gamma functions of the general form: Y(a 7 rf) = t a ~ x e l dt , [where 

Jjj 

a = \T iC , r/ = i(BDJ; + B)l D ], and is quite cumbersome. Because a very accurate alternative simplified 
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deflection equation could be introduced (see the following section) no attempt was made to program and use 
the cumbersome Lagrangianm deflection equation (49). 

Simplified Deflection Equation 

If A L t is replaced with A/, the simplified Eulerian deflection equation (40) can be converted into the 
following simplified Lagrangian deflection equation: 

y x = y-(tan0 M + tan<9 z ) + ; (/' = 1,2,3 ,....,«) (50) 

in which {tan 6 i X , tan 9 i } are to be calculated from the slope equation (45). 

Like the simplified Eulerian deflection equation (40), the simplified Lagrangian deflection equation (50) 
can also calculate extremely accurate deflections if accurate input slope data were used. Therefore, the need 
to program and use the extremely cumbersome equation (49) is avoided. 

STRUCTURE USED FOR SHAPE PREDICTION ACCURACY STUDIES 

For the shape-prediction accuracy studies of the newly formulated slope and deflection equations, a long 
aluminum tapered cantilever tubular beam (fig. 5) was used. This structure has geometries listed in table 1. 

As shown in figure 5, the strain-sensing stations [indicated with s i (i = 1,2,3,....,??)] are equally spaced 
along the bottom strain-sensing line, with strain-sensing stations, { £ 0’ £ n}’ respectively located at the fixed 
end (x = 0) and free end(x = /) . An upward point load of P = 300 lb (or P = 600 lb) was applied at the 
beam free end (fig. 5). 


Finite-Element Analysis 

For the shape prediction accuracy studies, reference slopes and deflections are needed as yardsticks. 
Therefore, finite-element analyses were first performed using SPAR (Structural Performance And Resizing) 
finite-element computer program (ref. 18) to generate the reference slopes and deflections. The SPAR 
program was developed by National Aeronautics and Space Administration Langley Research Center 
(Hampton, Virginia) in 1978 for preflight and subsequent flights (STS-l-STS-5) reentry heat transfer 
analyses of the Space Shuttle Orbiter Columbia until its loss on February 1, 2003 during STS- 107 reentry 
flight (12,500 miles per hour) at 203,000 feet altitude above North Central Texas. 

The structural part of SPAR can handle only linear analysis. The SPAR model was used to analytically 
generate surface strains for input to the Displacement Transfer Functions to calculate slopes and deflections 
for comparison with SPAR-calculated reference slopes and deflections for shape prediction accuracy 
analysis. Since aircraft wings are practically operating within linear range, the current technical publication 
is for linear analysis only. Structural buckling, collapsing failures, and geometrical nonlinearity were not 
considered. Keep in mind that small strains can produce large deflections for long span wings such as 
Global Observer (see: “Experimental Validations of Shape-Prediction Accuracies” section). 

Figure 6 shows the undeformed and deformed shapes of the SPAR model generated for the tapered 
cantilever tubular beam subjected to beam-tip load of P = 3001b (or P = 600 lb). In the figure, the size of 
the SPAR model is also indicated. Note that, for the loading case of P = 300 lb, the beam-tip deflection 
reached y n ~ 141.6 in, and the beam-tip slope reached as large as 0 n ~ 48 deg. For the doubled loading case 
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of P = 600 lb, the beam-tip deflection was doubled to y n ~ 283.2 in (because of linear elasticity), however, 
the beam-tip slope was not doubled, but only increased by 1.375 times to 0 n ~ 66 deg (= 1.375 X48 deg). 

Surface Strains 

For the present technical publication, the surface strains were analytically calculated by converting the 
nodal stresses or element stresses of the SPAR outputs into strains through stress-strain law. Because of 
finite-element modeling, the strain-sensing stations shown in figure 5 can be coincidental with the SPAR 
model nodal points or lying between the two adjacent nodal points along the bottom generatrix. 

Figure 7 shows the plots of SPAR-generated surface strains for the two loading cases P = {300, 600} lb. 
Note that the strains for the two loading cases increased practically linearly from the fixed ends and reached 
their respective peaks in the outboard region, then decreased rapidly toward the beam tip (because of 
decreasing depth factor), and finally became nearly zero at the beam tip. Theoretically, the strain at the 
beam tip (free end) should be zero. However, because of the finite-element discretization, the 
SPAR-generated beam-tip strain for each loading case showed a slightly non-zero value, which is 
proportional to the applied load, P, due to linear elasticity. Keep in mind that for the uniform cantilever 
beam, the strain curve is a straight line with maximum strain at the fixed end and zero strain at the free end 
(refs. 1-5). 

The SPAR-generated strain data of figure 7 were then input to the newly developed displacement 
transfer functions for the shape predictions. The results were then compared with the corresponding SPAR 
reference values for the shape prediction accuracy analysis. 

COMPARISONS OF SHAPE PREDICTION ACCURACIES 

In the comparisons of the shape prediction accuracies, the Shifted slope and deflection equations 
[eqs. (15)— (21)]; the Eulerian slope and deflection equations [eqs. (29)-(40)]; and the Lagrangian slope and 
deflection equations [eqs. (45)— (50)] were used. The slope and deflection prediction errors are respectively 
defined as the percent differences (deviations) from the SPAR-generated slope and deflection data, which 
were used as reference yardsticks. 


Slope Predictions 

Tables 2 and 3 respectively compare the slopes predicted from Shifted slope equations (15a)-(21a), the 
Eulerian slope equation (29), Lagrangian slope equation (45) for the loading cases of P = 300 lb and 
P = 600 lb. 
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Table 2. Comparisons of slopes calculated from SPAR and from different slope equations; long tapered tube 
( / = 300 in, c 0 =4 in, c n = 1 in) subjected to tip load ofP = 300 lb; n = 16; A / = 18.75 in. 


Slope, deg 

do 

e 2 

04 

06 

08 

010 

012 

014 

016 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

(% diff.) 

SPAR 

(reference) 

0.0000 

4.3526 

9.4255 

15.299 

21.9903 

29.3683 

37.0656 

44.2216 

48.2679 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 

(0.0000) 


Based on Shifted curvature equation (11) 


Non- 
uniform 
eq. (15a) 

0.0000 

(0.0000) 

4.3461 

(0.1493) 

9.4213 

(0.0466) 

15.2919 

(0.0523) 

21.9831 

(0.0327) 

29.3545 

(0.0470) 

37.0367 

(0.0780) 

44.1651 

(0.1278) 

48.0961 

(0.3559) 

Slightly 

0.0000 

4.3405 

9.4084 

15.2690 

21.9481 

29.3059 

36.9773 

44.1141 

48.1723 

tapered 
eq. (16a) 

(0.0000) 

(0.2780) 

(0.1814) 

(0.2020) 

(0.1919) 

(0.2125) 

(0.2382) 

(0.2431) 

(0.1981) 

1 st order eq. 
(17a) 

0.0000 

(0.0000) 

4.3426 

(0.2297) 

9.4131 

(0.1316) 

15.2771 

(0.1433) 

21.9594 

(0.1405) 

29.3191 

(0.1675) 

36.9864 

(0.2137) 

44.0967 

(0.2824) 

48.0164 

(0.5211) 

2 nd order 

0.0000 

4.3459 

9.4210 

15.2912 

21.9820 

29.3526 

37.0337 

44.1603 

48.0899 

eq. (18a) 

(0.0000) 

(0.1539) 

(0.0477) 

(0.0569) 

(0.0377) 

(0.0535) 

(0.0861) 

(0.1386) 

(0.3688) 

Improved 

0.0000 

4.3447 

9.4201 

15.2903 

21.9837 

29.3607 

37.0588 

44.2497 

48.2910 

eq. (19a) 

(0.0000) 

(0.1815) 

(0.0573) 

(0.0627) 

(0.0300) 

(0.0259) 

(0.0183) 

(0.0635) 

(0.0479) 

Log- 

expanded 
eq. (20a) 

0.0000 

4.3444 

9.4196 

15.2892 

21.9822 

29.3590 

37.0582 

44.2597 

48.3375 

(0.0000) 

(0.1884) 

(0.0626) 

(0.0699) 

(0.0368) 

(0.0317) 

(0.0200) 

(0.0862) 

(0.1442) 

Depth- 

0.0000 

4.3446 

9.4198 

15.2897 

21.9826 

29.3589 

37.0560 

44.2461 

48.2888 

expanded 
eq. (21a) 

(0.0000) 

(0.1838) 

(0.0605) 

(0.0667) 

(0.0350) 

(0.0320) 

(0.0259) 

(0.0544) 

(0.0433)* 


Based on Eulerian curvature equation (3) 


Eulerian 
non- 
uniform 
eq. (29) 

0.0000 

4.3544 

9.5072 

15.6656 

23.1294 

32.2247 

43.2335 

55.6503 

63.8446 

(0.0000) 

(0.0414) 

(0.8668) 

(2.3902) 

(5.1800) 

(9.7261) 

(16.6405) 

(25.8442) 

(32.2713) 

Lagrangian 
non- 
uniform 
eq. (45) 

0.0000 

4.3419 

9.3789 

15.1108 

21.4453 

28.0683 

34.4098 

39.5432 

41.9122 

(0.0000) 

(0.2458) 

(0.4944) 

(1.2360) 

(2.4784) 

(4.4265) 

(7.1651) 

(10.5794) 

(13.1679) A 
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Table 3. Comparisons of slopes calculated from SPAR and from different slope equations; long tapered tube 
( / = 300 in, c 0 =4 in, c n = 1 in) subjected to tip load ofP = 600 lb; n = 16; A / = 18.75 in. 


Slope, deg 

0« 

(% diff.) 

e 2 

(% diff.) 

04 

(% diff.) 

06 

(% diff.) 

0g 

(% diff.) 

010 

(% diff.) 

012 

(% diff.) 

014 

(% diff.) 

016 

(% diff.) 

SPAR 

(reference) 

0.0000 

(0.0000) 

8.6553 

(0.0000) 

18.3669 

(0.0000) 

28.6842 

(0.0000) 

38.9263 

(0.0000) 

48.3787 

(0.0000) 

56.4976 

(0.0000) 

62.8071 

(0.0000) 

65.9637 

(0.0000) 


Based on Shifted curvature equation (11) 


Non- 
uniform 
eq. (15a) 

0.0000 

(0.0000) 

8.6426 

(0.1467) 

18.3691 

(0.0120) 

28.6713 

(0.0450) 

38.9162 

(0.0259) 

48.3626 

(0.0333) 

56.4700 

(0.0489) 

62.7611 

(0.0732) 

65.8349 

(0.1953) 

Slightly 
tapered 
eq. (16a) 

0.0000 

(0.0000) 

8.6319 

(0.2704) 

18.3351 

(0.1731) 

28.6335 

(0.1768) 

38.8669 

(0.1526) 

48.3061 

(0.1501) 

56.4131 

(0.1496) 

62.7196 

(0.1393) 

65.8921 

(0.1085) 

1 st order 
eq. (17a) 

0.0000 

(0.0000) 

8.6360 

(0.2230 

18.3439 

(0.1252) 

28.6468 

(0.1304) 

38.8829 

(0.1115) 

48.3215 

(0.1182) 

56.4218 

(0.1342) 

62.7054 

(0.1619) 

65.7749 

(0.2862) 

2nd order 
eq. (18a) 

0.0000 

(0.0000) 

8.6424 

(0.1490) 

18.3585 

(0.0457) 

28.6703 

(0.0485) 

38.9147 

(0.0298) 

48.3605 

(0.0376) 

56.4672 

(0.0538) 

62.7572 

(0.0794) 

65.8303 

(0.2022) 

Improved 
eq. (19a) 

0.0000 

(0.0000) 

8.6400 

(0.1768) 

18.3569 

(0.0544) 

28.6686 

(0.0544) 

38.9171 

(0.0236) 

48.3698 

(0.0184) 

56.4911 

(0.0115) 

62.8300 

(0.0365) 

65.9811 

(0.0264) 

Log- 

expanded 
eq. (20a) 

0.0000 

(0.0000) 

8.6394 

(0.1837) 

18.3559 

(0.0599) 

28.6669 

(0.0603) 

38.9150 

(0.0290) 

48.3678 

(0.0225) 

56.4907 

(0.0122) 

62.8381 

(0.0494) 

66.0159 

(0.0791) 

Depth- 
expanded 
eq. (21a) 

0.0000 

(0.0000) 

8.6397 

(0.1802) 

18.3563 

(0.0577) 

28.6676 

(0.0579) 

38.9156 

(0.0275) 

48.3678 

(0.0225) 

56.4885 

(0.0161) 

62.8271 

(0.0318) 

65.9795 

(0.0240)* 


Based on Eulerian curvature equation (3) 


Eulerian 










non- 

0.0000 

8.7088 

19.0144 

31.3314 

46.2593 

64.4514 

86.4839 

86.2696 

88.1214 

uniform 
eq. (29) 

(0.0000) 

(0.6181) 

(3.5254) 

(9.2288) 

(18.8382) 

(33.2227) 

(53.0754) 

(37.3564) 

(33.5907) 


Based on Lagrangian curvature equation (7) 


Lagrangian 

non- 

0.0000 

8.6099 

18.0458 

27.4737 

35.8469 

42.0569 

44.9455 

44.9390 

44.9879 

uniform 
eq. (45) 

(0.0000) 

(0.2461) 

(1.7483) 

(4.2201) 

(7.9108) 

(13.0673) 

(20.4471) 

(28.4492) 

(31.7990) A 


< Acceptable range >|<- Poor range >| 

* Most accurate at beam tip A Worst case 


Note from tables 2 and 3 that the slope data predicted from the Shifted slope equations (15a) - (20a) are 
amazingly accurate beyond expectation. For the loading case of P = 300 lb (table 2), the beam-tip slope 
prediction errors are in the negligible range of 0.0433 to 0.5211 percent, and for the loading case of 
P = 600 lb (table 3), the beam-tip slope prediction errors decreased to the range of 0.0240 to 0.2862 percent 
(a nearly 45 percent reduction). The Depth-expanded slope equation (21a) produced the most accurate 
beam-tip slope data with infinitesimal prediction errors of only {0.0433, 0.0240} percent respectively for 
the loading cases of P = 300 lb (table 2) and P = 600 lb (table 3). 

In tables 2 and 3, the slopes calculated from the Eulerian slope equation (29) are referred to the 
deformed x-coordinate and are therefore somewhat larger than the corresponding SPAR slopes, which are 
referred to the un-deformed x-coordinate [that is, 6{s)> 6{x) , (fig. 2)]. The numbers listed in the 
parentheses under the slope data based on Eulerian formulation do not represent prediction errors, but are 
the percent differences because two different x-systems were used. When the load was doubled from 
P = 300 lb to P = 600 lb, the percent slope differences at the beam-tip based on Eulerian formulation 
increased from 32.2713 to 33.5907 percent (a nearly 4 percent increase). 
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The Lagrangian slope equation (45) predicted acceptable slope data in the inboard region, but the 
predicted slopes became lower than the corresponding SPAR reference slopes in the outboard region. When 
the load was doubled from P = 300 lb to P = 600 lb, the percent slope differences (errors) based on 
Lagrangian formulation increased from 13.1676 to 31.7990 percent (a nearly 141 percent increase) at the 
beam tip. 


Slope Curves 

The slope data listed in tables 2 and 3 are plotted in figures 8 and 9 respectively for the loading cases of 
P = {300, 600} lb for visual comparisons of the slope curves based on different formulations. As shown in 
figures 8 and 9, the slope curves based on the Shifted formulation [eqs. (15a) -(21a)] almost formed a single 
curve, which practically fell on top of the corresponding SPAR reference slope curve for the two loading 
cases P = {300, 600} lb, showing the high degrees of accuracies of the Shifted slope equations (15a)-(21a). 

In figure 8 (loading case, P = 300 lb), the slope curve of Eulerian formulation [eq. (29)] practically lies 
on top of the SPAR reference slope curve only in the inboard region, but branched out upward away from 
the SPAR reference slope curve in the outboard region because SPAR slopes are based on the undeformed 
v-system. Keep in mind that the Eulerian slopes [ 0(.y)] are always larger than the corresponding Lagrangian 
slopes [ 9(x) ] [that is, 9(s) > 9(x) , (fig. 2)]. 

The slope curve of Lagrangian formulation [eq. (45)] is quite close to the SPAR reference slope curve 
only in the inboard region, but branched out downward slightly away from the SPAR reference slope curve 
in the outboard region (fig. 8). 

In figure 9 (loading case, P = 600 lb), the slope curve based on the Eulerian formulation [eq. (29)], 
except for a small inboard region, continued to deviate upward further away from the SPAR reference slope 
curve in the outboard region. Keep in mind that the Eulerian slopes [#(*?)] are always larger than the 
corresponding Lagrangian slopes [ 9{x) ] [that is, 9{s) > 9{x) , (fig. 2)]. The peculiar wavy behavior of the 
Eulerian slope curve near the beam tip region could be attributed to the mathematical functional behaviors 
of equations (29) and (24) at large slopes approaching 90 deg (table 3) as a result of the magnified effect of 

the term (dy / dxf . 

At increasing load, P = 600 lb (fig. 9), the slope curve based on the Lagrangian formulation [eq. (45)], 
except for the small inboard region, deviated downward further away from the SPAR reference slope curve 
in the outboard region, causing the slope prediction errors to be greatly magnified. The rather poor slope 
predictions at increasing load (increasing deflections) could be attributed to the intensified effect of the term 

(dy / dx) 2 in the Lagrangian curvature equation (7). 

Deflection Predictions 

Tables 4 and 5 respectively compare the deflections calculated from the Shifted deflection equations 
(15b, c)-(21b, c), from the Eulerian deflection equation (40), and from Lagrangian deflection equation (50), 
for the loading cases of P = 300 lb and P = 600 lb. 
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Table 4. Comparisons of deflections calculated from SPAR and from different deflection equations; long 
tapered tube ( l = 300 in, c 0 =4 in, c n = 1 in) subjected to tip load of P= 300 lb; n = 16; A / = 18.75 in. 


Deflection, 

in. 

yo 

(% diff.) 

Y2 

(% diff.) 

Y4 

(% diff.) 

Y6 

(% diff.) 

ys 

(% diff.) 

Yio 

(% diff.) 

y 12 

(% diff.) 

yi4 

(% diff.) 

yi6 

(% diff.) 

SPAR 

(reference) 

0.0000 

(0.0000) 

1.4029 

(0.0000) 

5.8910 

(0.0000) 

14.0737 

(0.0000) 

26.6907 

(0.0000) 

44.7190 

(0.0000) 

69.3198 

(0.0000) 

101.7093 

(0.0000) 

141.6016 

(0.0000) 

Based on Shifted curvature equation (11) 

Non-uniform 
eq. (15b, c) 

0.0000 

(0.0000) 

1.3871 

(1.1262) 

5.8746 

(0.2784) 

14.0501 

(0.1667) 

26.6660 

(0.0925) 

44.6814 

(0.0841) 

69.2690 

(0.0733) 

101.6065 

(0.1011) 

141.3056 

(0.2090) 

Slightly 
tapered 
eq. (16b, c) 

0.0000 

(0.0000) 

1.3859 

(1.2118) 

5.8679 

(0.3921) 

14.0321 

(0.2956) 

26.6282 

(0.2342) 

44.6117 

(0.2399) 

69.1504 

(0.2444) 

101.4226 

(0.2819) 

141.1004 

(0.3540) 

1 st order 
eq. (17b, c) 

0.0000 

(0.0000) 

1.3863 

(1.1833) 

5.8704 

(0.3497) 

14.0387 

(0.2487) 

26.6416 

(0.1840) 

44.6348 

(0.1883) 

69.1847 

(0.1949) 

101.4585 

(0.2466) 

141.0562 

(0.3852) 

2nd order 
eq. (18b, c) 

0.0000 

(0.0000) 

1.3871 

(1.1262) 

5.8744 

(0.2818) 

14.0496 

(0.1712) 

26.6650 

(0.0963) 

44.6794 

(0.0886) 

69.2649 

(0.0792) 

101.5985 

(0.1089) 

141.2901 

(0.2200) 

Improved 
eq. (19b, c) 

0.0000 

(0.0000) 

1.3865 

(1.1690) 

5.8732 

(0.3022) 

14.0476 

(0.1855) 

26.6632 

(0.1030) 

44.6810 

(0.0850) 

69.2812 

(0.0557) 

101.6760 

(0.0327) 

141.5683 

(0.0235)* 

Log- 

expanded 
eq. (20b, c) 

0.0000 

(0.0000) 

1.3865 

(1.1690) 

5.8731 

(0.3039) 

14.0471 

(0.1890) 

26.6620 

(0.1075) 

44.6787 

(0.0901) 

69.2772 

(0.0615) 

101.6684 

(0.0402) 

141.5643 

(0.0263) 

Depth- 
expanded 
eq. (21b, c) 

0.0000 

(0.0000) 

1.3865 

(1.1690) 

5.8731 

(0.3039) 

14.0472 

(0.1883) 

26.6622 

(0.1068) 

44.6789 

(0.0897) 

69.2773 

(0.0613) 

101.6686 

(0.0400) 

141.5561 

(0.0321) 

Based on Eulerian curvature equation (3) 

Eulerian 
non-uniform 
eq. (40) 

0.0000 

(0.0000) 

1.3965 

(0.4562) 

5.8862 

(0.0815) 

14.0146 

(0.4199) 

26.4191 

(1.0176) 

43.7890 

(2.0797) 

66.7253 

(3.7428) 

95.3525 

(6.2500) 

128.0932 

(9.5397) 


Based on Lagrangian curvature equation (7) 


Lagrangian 
non-uniform 
eq. (50) 

0.0000 

1.3958 

5.8823 

14.0025 

26.3874 

43.7123 

66.5437 

94.9396 

127.5099 

(0.0000) 

(0.5061) 

(0.1477) 

(0.5059) 

(1.1364) 

(2.2512) 

(4.0048) 

(6.6559) 

(9.95 17) A 


< Acceptable range — >i< Poor range >| 

* Most accurate at beam tip A Worst case 
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Table 5. Comparisons of deflections calculated from SPAR and from different deflection equations; long 
tapered tube ( l = 300 in, c 0 =4 in, c n = 1 in) subjected to tip load ofP = 600 lb; n = 16; A / = 18.75 in. 


Deflection, 

in. 

yo 

(% diff.) 

Y2 

(% diff.) 

Y4 

(% diff.) 

Y6 

(% diff.) 

Ys 

(% diff.) 

Yio 

(% diff.) 

y 12 

(% diff.) 

y 14 

(% diff.) 

y i6 

(% diff.) 

SPAR 

(reference) 

0.0000 

(0.0000) 

2.8058 

(0.0000) 

11.7821 

(0.0000) 

28.1474 

(0.0000) 

53.3814 

(0.0000) 

89.4380 

(0.0000) 

138.6400 

(0.0000) 

203.4154 

(0.0000) 

283.2032 

(0.0000) 

Based on Lagrangian curvature equation (11) 

Log-Non- 
uniform 
eq. (15b, c) 

0.0000 

(0.0000) 

2.7742 

(1.1262) 

11.7492 

(0.2792) 

28.1001 

(0.1680) 

53.3319 

(0.0927) 

89.3628 

(0.0841) 

138.5379 

(0.0736) 

203.213 1( 
0.0995) 

282.6112 

(0.2090) 

Slightly 
tapered 
eq. (16b, c) 

0.0000 

(0.0000) 

2.7717 

(1.2153) 

11.7359 

(0.3921) 

28.0643 

(0.2952) 

53.2565 

(0.2340) 

89.2234 

(0.2399) 

138.3008 

(0.2447) 

202.8453 

(0.2803) 

282.2007 

(0.3540) 

1 st order 
Eq. (17b, c) 

0.0000 

(0.0000) 

2.7726 

(1.1833) 

11.7409 

(0.3497) 

28.0774 

(0.2487) 

53.2832 

(0.1840) 

89.2696 

(0.1883) 

138.3694 

(0.1951) 

202.9170 

(0.2450) 

282.1124 

(0.3852) 

2nd order 
eq. (18b, c) 

0.0000 

(0.0000) 

2.7742 

(1.1262) 

11.7489 

(0.2818) 

28.0993 

(0.1709) 

53.3300 

(0.0963) 

89.3587 

(0.0887) 

138.5299 

(0.0794) 

203.1971 

(0.1073) 

282.5803 

(0.2199) 

Improved 
eq. (19b, c) 

0.0000 

(0.0000) 

2.7731 

(1.1654) 

11.7465 

(0.3022) 

28.0952 

(0.1855) 

53.3264 

(0.1030) 

89.3619 

(0.0851) 

138.5624 

(0.0560) 

203.3521 

(0.0311) 

283.1367 

(0.0235)* 

Log- 

expanded 
eq. (20b, c) 

0.0000 

(0.0000) 

2.7731 

(1.1654) 

11.7462 

(0.3047) 

28.0943 

(0.1886) 

53.3240( 

0.1075) 

89.3573 

(0.0902) 

138.5543 

(0.0617) 

203.3368( 

0.0386) 

283.1285 

(0.0264) 

Depth- 
expanded 
eq. (21b, c) 

0.0000 

(0.0000) 

2.7730 

(1.1690) 

11.7462 

(0.3047) 

28.0944 

(0.1883) 

53.3244 

(0.1068) 

89.3579 

(0.0896) 

138.5545 

(0.0617) 

203.3372 

(0.0384) 

283.1123 

(0.0321) 

Based on Eulerian curvature equation (3) 

Eulerian 
non-uniform 
eq. (40) 

0.0000 

(0.0000) 

2.79213 

(0.4887) 

11.7129 

(0.5866) 

27.6060 

(1.9232) 

51.1125 

(4.2505) 

82.4108 

(7.8571) 

131.4861 

5.1598) 

130.7083 

(36.7442) 

151.6729 
(46.443 8) A 


Based on Lagrangian curvature equation (7) 


Lagrangian 

non-uniform 


eq. (50) 


0.0000 

( 0 . 0000 ) 


2.7865 

(0.6879) 


11.6815 

(0.8538) 


27.5015 

(2.2947) 


50.8036 

(4.8290) 


81.3751 

(9.0151) 


117.2992 

(15.3930) 


154.5417 

(24.0265) 


191.8715 

(32.2495) 


< Acceptable range >|<- 


■Poor range- 


>i 


* Most accurate at beam tip A Worst case 


Note from tables 4 and 5 that the deflections calculated from all the Shifted deflection equations 
(15b, c)-(21b, c) are amazingly accurate beyond expectations with negligible beam-tip prediction errors in 
the range of 0.0235 to 0.3852 percent for both loading cases, P = {300, 600} lb. Note also from tables 4 and 
5 that under linear analysis, the predicted deflections and associated predicted deflection errors (in inch) are 
linear functions of the applied load. Because of linearity, the deflection prediction errors (in percentage) are 
invariant to the loading magnitude. The improved deflection equations (19b, c) turned out to be the most 
accurate deflection equation, which produced only infinitesimal amount of prediction error of 
0.0235 percent at the beam-tip for both loading cases, P = {300, 600} lb. 


The Eulerian deflection equation (40) predicted acceptable deflection data in the inboard region, but the 
prediction errors continued to magnify toward the beam tip. Because the slope and deflections are not 
directly proportional to the surface strains in the Eulerian formulation [eqs. (34) and (40)], when the applied 
load, P = 300 lb, was doubled to P = 600 lb (surface strains doubled), the deflection prediction errors were 
greatly magnified from 9.5397 to 46.4438 percent (a nearly 387 percent increase) at the beam tip. 


The Lagrangian deflection equation (50) also predicted acceptable deflection data in the inboard region, 
beyond which the prediction errors continue to magnify toward the beam tip. Because the slope and 
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deflections are not directly proportional to the surface strains in the Lagrangian formulation [eqs. (45) and 
(50)], when the applied load, P = 300 lb, was doubled to P = 600 lb, the deflection prediction errors 
increased nonlinearly from 9.9517 to 32.2495 percent (a nearly 224 percent increase) at the beam tip. 

Deflection Curves 

The deflection data listed in tables 4 and 5 are plotted in figures 10 and 11 respectively for the loading 
cases of P = {300, 600} lb for visual comparisons of the deflection curves based on different formulations. 

As shown in figure 10 (loading case of P=300 lb), all the deflection curves based on the Shifted 
formulation [eqs. (15b, c)-(21b, c)] graphically fell on top of the associated SPAR reference deflection 
curves, showing exceptional accuracy of the Shifted deflection equations (15b, c)-(21b, c). 

Note from figure 10 that for loading case of P = 300 lb, the deflection curves based on the Eulerian 
formulation [eq. (40)] and Lagrangian formulation [eq. (50)] are practically coincidental, and agreed with 
the SPAR reference deflection curve only in the inboard region, but deviated slightly downward from the 
SPAR reference deflection curve in the outboard region (fig. 10). 

Figure 1 1 shows that at double load, P = 600 lb, the deflection curves based on the Shifted formulation 
[eqs. (15-b, c)-(21-b, c)] are again extremely close to the SPAR-reference deflection curve, showing 
amazing high degrees of prediction accuracies (tables 4, 5). The deflection curve of the Eulerian formulation 
agreed with the SPAR reference deflection curve only in the inboard region and deviated downward from 
the SPAR reference deflection curve in the outboard region, and became wavy near the beam tip, showing 
poor prediction performances. Again, the peculiar wavy behavior of the Eulerian deflection curve near the 
beam tip region is associated with the wavy behavior of the slope curve shown in figure 9, and could be 
attributed to the mathematical functional behaviors of equations (40) and (24) at large slopes approaching 
90 deg (table 3). 

As shown in figure 11, at doubled load P = 600 lb, the deflection curves based on the Lagrangian 
formulation [eq. (50)] fell on top of the SPAR reference deflection curve only in a small inboard region, and 
then branched out downward further away from the SPAR reference deflection curve toward the beam 
outboard region, giving discouraging shape predictions. 

SHAPE PREDICTION ERROR CURVES 

The slope and deflection predictions errors (percent differences from associated SPAR data) listed in the 
parentheses in tables 2-5 are plotted in the following sections for visual comparisons. 

Slope-Error Curves 

Figure 12 shows the slope-error curves for P = 300 lb loading case associated with different 
Displacement Transfer Functions. The slope differential curve based on Eulerian formulation shown in 
figure 12 is not the slope error curve, but the curve showing the difference in slopes referred to the deformed 
and un-deformed v-coordinate. The slope differential for the Eulerian formulation steadily increased 
span-wisely and reached as large as 32.2713 percent at the beam tip (table 2). 

The slope-prediction errors based on the Lagrangian formulation continued to increase almost linearly 
span-wisely and reached a peak value of 13.1676 percent at the beam tip (table 2). 
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The slope-error curves (fig. 12) based on the Shifted formulation practically formed a single curve along 
the 0 percent-error line and was indistinguishable; therefore, those slope error curves had to be re -plotted 
independently with a different scale. 

Figure 13 shows the slope-error curves for the P = 300 lb loading case, calculated from seven Shifted 
slope equations (15a)-(21a). For the present example structure, the Depth-expanded case [eq. (21a)] 
produced the lowest slope-prediction errors near the beam tip region, with a slope prediction error of merely 
0.0433 percent at the beam tip (table 2). 

Figure 14 shows the slope-error curves for the P = 600 lb loading case associated with nine different 
displacement transfer functions. The lofty curve of the Eulerian formulation shown in figure 14 is not the 
slope-error curve, but a slope-differential curve for showing the difference in slopes based on the deformed 
x-coordinate (Eulerian) and un-deformed x-coordinate (SPAR). The peak slope differential of 
53.0754 percent occurred at the strain-sensing station i =12 (table 3). 

The slope-error curve of the Lagrangian formulation (fig. 12) is slightly ^-shaped and tilted away from 
the 0 percent-error axis with peak error of 31.7990 percent at the beam tip (table 3). The seven slope-error 
curves based on the Shifted formulation (fig. 14) are crowded together near the 0 percent-error axis and 
practically indistinguishable. Therefore, for better visual comparisons, those seven slope-error curves based 
on the seven Shifted formulations had to be re-plotted with a different scale. 

Figure 15 shows the slope-error curves for the P = 600 lb loading case calculated from the seven Shifted 
slope equations (15a)-(21a). For the present example structure, the Depth-expanded case [eq. (21a)] 
produced the lowest slope-error of 0.0240 percent at the beam tip ( table 3). Because the beam slopes are not 
directly proportional to the applied load, increasing the load from P = 300 lb to P = 600 lb caused the 
overall beam tip slope errors to come down by nearly 50 percent (compare figures 13 and 15). 

Deflection-Error Curves 

Figure 16 shows the deflection-error curves for the P = 300 lb loading case associated with nine 
different displacement transfer functions. Note that both the Eulerian and Lagrangian formulations produced 
similar concave upward deflection-error curves lying far above the 0 percent-error axis, with respective 
deflection prediction errors as large as {9.5397, 9.9517}percent at the beam tip (table 4). Because the 
deflection-error curves based on the Shifted formulation are so close to the 0 percent-error axis, and are 
almost indistinguishable in figure 16, those deflection-error curves were re -plotted with a different scale for 
better visual comparisons. 

Figure 17 shows the deflection-error curves for the P = 300 lb loading case, calculated from the seven 
sets of the Shifted deflection equations. (15b, c)-(21b, c)]. For the present example structure, the Improved 
case [eq. (19b, c)] produced the lowest deflection-prediction errors in the outboard region, with negligible 
deflection-prediction error of 0.0235 percent at the beam tip (table 4). The increasing deflection errors at 
strain-sensing station, i = 2, is due to a small predicted deflection error divided by a small SPAR reference 
deflection. 

Figure 18 shows the deflection-error curves for the P = 600 lb loading case associated with nine 
different displacement transfer functions. The deflection errors based on the Shifted formulation are very 
small as compared to those of both the Eulerian and Lagrangian formulations. The dent in the Eulerian 
deflection-error curve (fig. 18) is associated with the hump in the Eulerian slope-error curve (fig. 12). 
Because the deflection-error curves based on the Shifted formulation are so close and are almost 
indistinguishable in figure 18, those deflection-error curves were re -plotted using an enlarged scale for 
better visual comparisons. 
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Figure 19 shows the deflection-error curves for the P = 600 lb loading case calculated from the seven 
sets of Shifted deflection equations (15b, c)-(21b, c). Because the beam deflections are directly proportional 
to the strains (or applied load), by increasing the load from P = 300 lb to P = 600 lb, the overall deflection 
errors remain practically the same (tables 4 and 5). Again, the increasing deflection errors at strain-sensing 
station, i = 2 , is due to a small deflection error divided by a small SPAR reference deflection. For the 
present example structure, the Improved case [eq. (19b, c)] produced the lowest deflection errors in the 
outboard region, with a negligible deflection prediction error of 0.0235 percent at the beam tip (same as the 
P = 300 lb case) (that is, figs. 17 and 19 are practically identical). 

EXPERIMENTAL VALIDATIONS OF SHAPE-PREDICTION 

ACCURACIES 

In the process of writing the present technical publication, large-scale ground loads tests were conducted 
on a long span (high aspect ratio) full-scale tapered wing (175-ft wingspan) at the NASA Dryden Flight 
Loads Laboratory. These ground tests were used to validate the shape prediction accuracies of the 
Displacement Transfer Functions. The results were reported in a separate NASA technical publication 
authored by Jutte, Ko, et al (ref. 19). 

In the ground tests, the whole wing was loaded from 0 percent up to 100 percent design limit load 
(DLL) at which the wing-tip slope exceeded 20 deg. Surfaces strains were measured by using four 
fiber-optic strain-sensing lines installed on the wing surfaces (two on the wing upper surface, two on the 
wing lower surface). Wing deflections were measured by means of a photogrammetry system, and were 
used as yardsticks to evaluate the accuracies of the predicted wing deflections. Using the measured surface 
strains as inputs, theoretical wing deflections were calculated from the recursive Displacement Transfer 
Functions [eqs. (16a, b)], which were formulated for slightly tapered beams (based on piecewise-linear 
strain representations) such as the current slightly tapered wing structure. 

For completeness of the present technical report, wing deflection plots in figures 11, and 12 of reference 
19 were duplicated respectively as figures 20, and 21 in the present technical report. Figures 20, and 21 
respectively show the plots of the predicted and experimentally measured wing deflections along the 
number 3 (FWD) and number 4 (AFT) strain-sensing lines. Amazingly, the measured data points practically 
fell on top of the associated predicted deflection curves for the whole range of loading levels 
[(0-100) percent DLL]. The wing-tip prediction errors associated with the number 3 and number 4 
strain-sensing lines are respectively in the negligible ranges of (0.09-1.72) percent and (0.25-1.45) percent 
at 100 percent DLL depending on the sensor-to-sensor distances used. For details of the experimental 
shape-prediction studies, please see reference 19. 

Note from figures 20 and 21 that the maximum wing-tip deflections of the Global Observer reached 
approximately 160 in. Hence, the ratio of (maximum wing-tip deflection) to (half wingspan) can be 
calculated as: 


Maximum wingtip deflection 
Half wingspan 


160 

(175/2) x 12 


= 0.1524 = 15.24% 


(51) 


The 15.24 percent normalized wing-tip deflection from equation (51) indicates that the ground test of 
the Global Observer was operating in a borderline region between the linear and nonlinear regimes. 

Keep in mind that each Shifted Displacement Transfer Function was formulated by geometrically 
relating the surface strains to curvatures, and then to deflections. Therefore, one could expect that the 
Shifted Displacement Transfer Function can faithfully compute the right deflections whether the input 
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surface strains come from linear or nonlinear analysis. To find out whether a good correlation can be 
obtained with realistic truly large deformation cases, nonlinear analysis will be carried out separately, and 
the results will be published in a subsequent NASA publication. 

The Global Observer ground test results (ref. 19) thus validated the mathematical accuracies of the 
Shifted Displacement Transfer Functions, and established confidence in using the Shifted 
Displacement Transfer Functions for large-deformation shape predictions. 

DISCUSSIONS 

The vitally important lesson learned in the present technical publication is that for large deformations, 
one cannot just naively integrate the classical Eulerian and Lagrangian curvature equations {(12), (13)} 
(written in the x-y system) without considering the true nature of large deformations. After spending long 
hours in the lengthy mathematical derivations because of keeping the nonlinear term (< dy/dx ) 2 , both Eulerian 
and Lagrangian formulations resulted in poor shape predictions at large deformations, and contrary to the 
hopeful expectation, no better shape prediction accuracies could be achieved. 

In the piecewise formulations, the embedded beam is discretized into multiple small domains. Within 
each small domain, the deformations are relative to the inboard domain juncture, and can be considered 
small. Therefore, in the piecewise integrations of Eulerian or Lagrangian curvature equations {(12), (13)}, 
the nonlinear term (< dy/dx ) 2 , could be neglected to achieve high accuracy shape predictions like the Shifted 
case. 

For large bending of the beam (for example, fishing pole), the true deflection, y , is actually the curved 
distance traced by the same material point from its initial undeformed position to its final deformed position 
(fig. 3). Thus, the conventional vertical deflection, y, is merely the vertical component of the curved 
deflection, y . For small deformations, the curved deflection, y , can be approximated with the vertical 
deflection, y . However, for large deformations, the difference between y and y becomes larger (y > y). 
For large deformations, one could write the curvature equation in the s- y system (instead of x-y system) as 
equation (52): 


R(s ) ds ds\ds) ds 2 ds 


(52) 


which contains no classical nonlinear term (dy/dx) 2 (see also Appendix B). 

In view of the incremental length relationship, ds 2 = dx 2 +dy 2 (fig. 1), the original Eulerian curvature 
equation (3) can be rewritten in the following revised forms for the in-extensional beam (ds = dx) as: 


a. Revised Eulerian Curvature equation (changing deformed x-coordinate into undeformed x-coordinate, 
and then setting ds = dx): 


R 


d 2 y/dx 2 d 2 y {dx) 2 

[1 +{dyldxff 2 ~ dx 2 (ds? 

\ v / 

In deformed x -system 


d 2 y dx 
ds 2 ds 


Deformed x 
undeformed x, 
then set dx=ds 


dx 2 ' 

Shifted 


(53) 
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Also, in view of the incremental length relationship, ds 2 = ( dx-du ) 2 + dy 2 (fig. 2), the original 
Lagrangian curvature equation (4) can be rewritten in the following revised forms for the in-extensional 
beam (ds = dx) as: 

b. Revised Lagrangian curvature equation (setting du = 0,ds = dx) 

1 d 2 y/d.s 2 d 2 y ds cfy ds d 2 y ds 

R )l-(dyidsf ^ )ds 2 -dy 2 & 2 ) (dx-du) 2 <& 2 &~du 

' v ' 

In undefonned x -system 

•Py ds ( 54) 

ds 2 d(x-u) 

Shifted 

Note that, for the in-extensional beams (ds = dx ), both equations (53) and (54), written in undeformed 
x-system (x = s ), degenerated into the Shifted curvature equation (11), containing no nonlinear term 
(dy / dx) 2 . Thus, the formulations based on the revised Eulerian and Lagrangian curvature equations 
{(53) or (54)}, become the Shifted formulation, and can provide accurate shape predictions. 

CONCLUDING REMARKS 

Large-Deformation Displacement Transfer Lunctions were formulated through piecewise integrations of 
the Eulerian curvature-strain equation, the Lagrangian curvature-strain equation, and the Shifted 
curvature-strain equation. A long tapered cantilever tubular beam was chosen in the shape-prediction 
accuracy studies of the newly developed Displacement Transfer Lunctions. Surface strains for input to the 
Displacement Transfer Lunctions were analytically calculated from the finite-element analysis of the tapered 
cantilever tubular beam. Also, the finite-element-generated slopes and deflections were used as reference 
yardsticks in the prediction accuracy studies. The results are itemized below: 

1 . Changing from piecewise-linear strain representations to piecewise-nonlinear strain representations, 
the shape prediction accuracies could be improved considerably for the tapered cantilever tubular beam. 

2. Based on the Shifted curvature equation, the resulting seven sets of the Shifted Displacement Transfer 
Functions were found to be amazingly accurate beyond expectation for large-deformation shape predictions, 
even for the case with a beam-tip slope reaching as large as 66 deg. The beam-tip shape prediction errors of 
the seven sets of the Shifted slope and deflection equations, are in the negligible ranges of 
{0.0240 - 0.521 1 } percent for the slopes, and {0.0235 - 0.3852} percent for the deflections. 

3. For the Shifted formulation (based on linear analysis), the deflections are directly proportional to the 
surface strains, and therefore, for the same type of loading condition, the negligible deflection prediction 
errors are theoretically independent of the loading magnitude (that is, degree of deformation). 

4. The amazing accuracy of the Shifted formulation could be attributed to: 1) shifting of the deformed 
material points back to their respective initial undeformed x-positions; and 2) piecewise integrations of the 
Shifted curvature equation, enabled by dicretization of the embedded beam. 

5. Based on the classical Eulerarian curvature equation, the resulting slope and deflection equations 
gave poor shape predictions at large deformations because of retaining the term (dy / dx) 2 in the discretized 
approach. The beam tip deflection prediction errors are in the poor ranges of {9.5397 - 46.4438} percent. 
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6. Based on Eulerian formulation, when the applied load was doubled (that is, doubling deflections), 
because the deflections are not directly proportional to the surface strains, the deflection prediction errors 
were greatly magnified from 9.5397 percent to 46.4438 percent at the beam-tip (a 386.8476 percent 
increase). 

7. Based on the physical (Lagrangian) curvature equation, the resulting slope and deflection equations 
gave discouraging shape predictions at large deformations because of retaining the term (< dy / dx) 2 in the 
discretized approach. The beam-tip prediction errors are in the poor range of {13.1676 - 31.7990} percent 
for the slopes, and {9.9517 - 32.2495} percent for the deflections. 

8. Based on the Lagrangian formulation, when the applied load was doubled (that is, doubling 
deflections), because the slopes and deflections are not directly proportional to the surface strains, the 
beam-tip slope prediction errors were magnified from 13.1676 percent to 31.7990 percent (a 141.4943 
percent increase), and the beam-tip deflection prediction errors increased from 9.9517 percent to 32.2485 
percent (a 224.0602 percent increase). 

9. If the deformed material points of an in-extensional beam are shifted back to their respective initial 
undeformed x-locations, then the term {dy / dx) 2 in the Lagrangian curvature equation will become zero, 
implying linearization. The resulting linearlized Lagrangian formulation will become the Shifted 
formulation, and provides equally accurate shape predictions. 

10. The linearized Eulerian curvature equation [the term {dy / dx) 2 neglected], even with similar 
mathematical form, cannot be equivalent to the Shifted curvature equation because, the former is referred to 
the deformed x-coordinate, while the latter to the undeformed x-coordinate. The two x-systems will become 
quite different at large deformations. 

11. The simplified alternative deflection equations developed for the Eulerian and Lagrangian 
formulations could not provide good shape predictions at large deformations because of the poor input slope 
data at large deformations. 

12. Retaining the term {dy / dx) 2 in the piecewise integrations of curvature equations resulted in rather 
poor shape prediction performances of both the Eulerian and Lagrangian formulations. 

13. For in-extensional beams {ds = dx ), by writing the original Eulerian and Lagrangian curvature 
equations in 5-system, the formulations based on the revised Eulerian and Lagrangian curvature equations 
can degenerate into the Shifted formulation and, therefore, could provide accurate shape predictions. 

14. The Global Observer ground test results validated the mathematical accuracies of the Shifted 
Displacement Transfer Functions formulated, and giving confidence in using the Displacement Transfer 
Functions for large-deformation shape predictions. 
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FIGURES 
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Figure 1. Small segment of the deformed nonuniform embedded beam for geometrically relating local radius 
of curvature, R(s ), to associated surface bending strain. 



Figure 2. Small segment of the deformed-beam elastic curve used for deriving the physical (Lagrangian) 
curvature equation. 
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Figure 3. Small segment of the deformed-beam elastic curve used for deriving the Shifted curvature 
equation. 



Figure 4. Geometry used for calculations of the projected domain length A L i? along with /-axis in the 
Eulerian description. 
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Figure 5. Tapered cantilever tubular beam with strain-sensing stations evenly distributed along the bottom 
strain-sensing line. 


Size of SPAR model 



Figure 6. Large deformation of the tapered cantilever tubular beam (SPAR finite-element model) subjected 
to a vertical tip load of P = 300 lb (or P = 600 lb); / = 300 in; c 0 = 4 in; c„ = 1 in. 
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Figure 7. SPAR-generated surface-strain curves for the tapered cantilever tubular beam subjected to tip 
loads of P = {300,600} lb; / = 300 in; c 0 = 1 in. 



Figure 8. Slopes of the tapered cantilever tubular beam calculated from slope equations based on the 
Eulerian, Lagrangian, and Shifted formulations; / = 300 in; c 0 = 4 in; c„ = 1 in; P = 300 lb. 
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Figure 9. Slopes of the tapered cantilever tubular beam calculated from slope equations based on the 
Eularian, Lagrangian, and Shifted formulations; / = 300 in; c 0 = 4 in; c„ = 1 in; P = 600 lb. 



Figure 10. Deflections of the tapered cantilever tubular beam calculated from deflection equations based on 
the Eularian, Lagrangian, and Shifted formulations; / = 300 in; c 0 =4 in; c„ = 1 in; P = 300 lb. 
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Figure 11. Deflections of the tapered cantilever tubular beam calculated from deflection equations based on 
the Eulerian, Lagrangian, and Shifted formulations; / = 300 in; c 0 =4 in; c„= 1 in; P = 600 lb. 
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Figure 12. Comparisons of 9j prediction errors of different displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 300 lb at tip. 
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Figure 13. Comparisons of 6 t prediction errors of Shifted displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 300 lb at tip. 
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Figure 14. Comparisons of Q x prediction errors of different displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 600 lb at tip. 
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Figure 15. Comparisons of 6 t prediction errors of Shifted displacement transfer functions; long tapered 
cantilever tubular beam, / = 300 in; c 0 = 4 in; c n = 1 in; P = 600 lb at tip. 



Figure 16. Comparisons of y x prediction errors of different displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 300 lb at tip. 
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P = 300 lb 


Figure 17. Comparisons of y t predicition errors of Shifted displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 300 lb at tip. 
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Figure 18. Comparisons of y t prediction errors of different displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 600 lb at tip. 
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P = 600 lb 



Figure 19. Comparisons of y t prediction errors of Shifted displacement transfer functions; long tapered 
cantilever tubular beam; / = 300 in; c 0 = 4 in; c n = 1 in; P = 600 lb at tip. 
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Figure 20. Comparison of predicted and measured wing deflections along the number 3 strain-sensing line 
for different loading levels (duplication of figure 11, ref. 18). 
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Figure 21. Comparison of predicted and measured wing deflections along the number 4 strain-sensing line 
for different loading levels (duplication of figure 12, ref. 18). 
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APPENDIX A 


ALTERNATIVE DERIVATION OF THE PHYSICAL (L AGRAN GIAN) 

CURVATURE EQUATION 

The alternative method used by Kopmaz, et al (ref. 12) and Hodges (ref. 13) in deriving the physical 
(Lagrangian) curvature equation (7) is described in the following paragraphs. 

In view of figure 2, when the beam is bent, the undeformed material points, {^4(x,0) , B(x + dxfi)} , will 
be moved to respective deformed positions, {A'[{x - u ), y] , B'[(x + dx ) ~(u + du ), (y + dy)]} , in which { u , 
y} respectively denote the horizontal and vertical displacement of material point A (moving to point A'). 
Then, in view of the curved triangle A'B'C , the infinitesimal curved length, ds{= A'B'), along the elastic 
curve of the deformed beam can be expressed in terms of {dx, du, dy } as: 


ds = ( dx - du) 2 + {dy) 2 

Dividing through equation (Al) by dx , one obtains equation (A2): 


dx V v dx 


du^\ ( dy^ 


+ 


\dx j 


(Al) 


(A2) 


In view of curved triangle A'B'C (fig. 2), the slope, tan 6{s) , at point A' can be expressed as equation 
(A3): 


tan 6{s) = — — — 
dx - du 


dy I dx 
1 -du / dx 


from which, one obtains equation (A4): 


9(s ) = tan 1 


r dy I dx ^ 
yl~ du/dXy 


Differentiating equation (A-4) with respect to x, one obtains equation (A5): 


d0(s) 

dx 


(d 2 y / dx 2 )( 1 - du/dx) + (dy/dx)(d 2 u/dx 2 ) 
(1 - du/dx) 2 

1 f dy/dx V 
yl - du/dx j 


d 2 y 


dx 2 


j du 
dx 


+ 


dy d 2 u 
dx dx 2 


j dll' 

2 

+ 


d y J 


ydXy 


The curvature, l/R , of the elastic curve of the deformed beam can be written as: 

1 _ dd(s) _ d9(s) dx 
R ds dx ds 


(A3) 


(A4) 


(A5) 


(A6) 
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Substitutions of equations (A2) and (A5) into equation (A6) yields equation (A7): 


d 2 y 

f 

1- 

du ^ 

dy 

+ — 

d 2 u 


d 2 y 

( 

1 

1 _ dx 2 


dx y 

dx 

dx 2 

i 

dx 2 

V 


du 

dx 


+ 


dy d 2 u 
dx dx 2 


R 


f du 
dx 


+ 


J 

~deTdT 


r dy } 2 
\dx j 



1- 


du 
dx . 


r 


+ 


dy 


v 


A 
dx ) 


) “ 1 3/2 


dx! ds 


(A7) 


For the in-extensional beam, the curved coordinate, s, of any deformed material point is always equal to 
its initial un-deformed x-coordinate (fig. 2). Namely, s = x and ds = dx , and equation (A2) can be written 
as: 


r ' du ^ 

1 

+ 

(dy\ 

V dX y 


K dX y 


Rearranging equation (A8) as equation (A9): 


r du V 
y dx y 


- du ( Jv Y 
2 — + — 
dx \dx J 


= 0 


Solving for duldx one obtains: 


— =i± \\-(— 
dx \ v dx y 


which can be rewritten as: 



dx 



>y 

K dXy 


Differentiating equation (A10), du 2 / dx 2 can be expressed as: 



dy dy 

_ J dx ofc 2 


r dy } 2 

\dx j 


(A8) 


(A9) 


(A10) 


(All) 


(A 12) 
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Substitutions of equations (All) and (A 12) into the curvature equation (A7) yields equation (A 13): 


~R 


cTy 

dx 2 


+ J1- 


dv 

dx 


+ 


dy d 2 y 
dy dx dx 2 


dx 


L 

i dy ) 

2 _ d 2 y 

+ — j- 

i- 

( dy ) 

2 

+ 

'dy \ 

i 

ydx j 

dx 


K dX y 



ydXy 


dx 2 


2 


+ Jl ' 


dy_ 
\dx j 


+ 


dy_ 

dx 


3/2 


1- 


dy ' I , f dy 


\dx j 


+ 


\dx j 


3/2 


l- 


\dx j 


(A13) 


After cancelling terms, and choosing the positive sign, one obtains equation (A 14): 

1 _ d 2 y / dx 2 

R sfl - {dyldxf (A14) 

which is the physical (Lagrangian) curvature equation (7) for the in-extensional beam. 

Keep in mind that dy / dx in equation (A 14) is referred to the undeformed x-system. From right triangle, 
A (A'DE) (fig. 2), the slope, d(x) , in the x-system can be expressed as equation (A15): 


#(x) = tan 1 


dy A 
dx j 


(A 15) 


Comparing equations (A 15) and (A4), one notes that #(x) in the undeformed x-system is slightly 
smaller than 6>(.y)[=sm x (dyl ds)\ in the 5-system [that is, 0{x) < 0(.s)] (fig. 2). 


In the bending deformations, the axial displacement, u , of a material point is very small. If the axial 
displacement, u , is neglected (that is, u —> 0 ) in equation (A4), then 0(s) will become 6{x) of equation 
(A15). The condition (u —> 0) implies that the deformed material points {A f ,5'} (fig. 3) are horizontally 
shifted back to points {A" , B "} at their respective initial undeformed x-locations. It is important to mention 
that u — » 0 will make dy/dx^O according to equation (A8). 
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APPENDIX B 


GRAPHICAL DERIVATION OF THE SHIFTED CURVATURE EQUATION 

For the in-extensional beam, the curved coordinate, a, of any deformed material point is always equal to 
its initial un-deformed coordinate x (fig. 3). Namely, s = X and ds — dx, then the basic curvature equation 
can be written as equation (Bl): 


1 _ dQ(s) _ dQ(s) 

R{s) ds dx 


(Bl) 


In the Eulerian or Lagrangian curvature equation (3) or (7), the term (dy I dx) 2 appeared because, the 
slope, 6{s ) , of the beam elastic curve was expressed in term of either \2mC l {dyldx) or sin ~\dy/ds ) . It 
was found that term (dy I dx) 2 caused the formulated Displacement Transfer Functions to give rather poor 
shape predictions at large deformations. To avoid the appearance of the term (dy / dx ) 2 , one can 
alternatively express the slope, 6{s ) , in term of the curved deflection, y , (fig. 3), which is the curved 
length AA' actually traced by the same material point from its initial un-deformed position A to its final 
deformed position A f . Thus, in view of the infinitesimal curved triangle, A'B'C' [ C' will converge to 
BB' -curve in the limit, d6 ^ 0 , (fig. 3)], the slope, 0(s ) , can then be expressed as: 




(B2) 


By combining equations (Bl) and (B2), one obtains the following curvature equation (B3) associated 
with the curved deflection, y , as: 


1 _ d 2 y _ d 2 y 

R(s) ds 2 dx 2 


(B3) 


In the formulations of the Shifted displacement transfer functions, the deformed material points 
{A f B'} (fig. 3) are shifted horizontally to points {A" ,B") with x-coordinates respectively matching the 
original un-deformed x-coordinates of points {A,B} . The curve connecting the shifted deformed material 
points will be slightly detached from the true elastic curve of the actual deformed beam. Keep in mind that, 
in the bending deformations, the axial displacements are miniscule, and so also the amounts of horizontal 
shifting. 

Graphically, the horizontal shifting (fig. 3) will cause the curved deflection, y, to distort into a straight 
vertical deflection, y, and also causing the curved triangle, A'B'C' , to distort into a right triangle, tsA!'B"C " . 
Thus, the shifted-slope, #(x)[=tan ~ { (dy/dx)\, in the x-system will be slightly smaller than the original 

slope, 0(s)[= sin ~ { (dy/ds)] in the A-system [that is, 0(x) < 6{s ) , (fig. 3)]. After the horizontal shifting, and 
ignoring the resulting small perturbation of R(x ), equation (B-3) can be rewritten as equation (B-4): 


1 _d 2 y 

RM dx* 


(B4) 
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which is similar to equation (11) in the main text, and is called the Shifted curvature equation for the 
in-extensional beam (ds = dx). 
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APPENDIX C 


DERIVATION OF PROJECTED DOMAIN LENGTH FOR THE 
EULERIAL FORMULATION 

In view of curvature-strain equation (2), the averaged curvature, 1 / R (fig. 4), for the elastic curve 
bounded by Ad can be approximated by using the averaged strains and the averaged depth 
factor as equation (Cl): 


_ (g,-! + g ,) /2 _ g, i +g, 

R (c m + C< )/2 c M+Ci 1 ’ 

In view of figure 4, for the in-extensional beam, one can relate the chord length, A l c , between points 
{ i - 1 , /'} , to the curved domain length, A / , as: 


A l c = 2 R sin 


A 0 
2 


= 2 R sin 


A l_ 
2 R 


(C2) 


The deformed domain length, AT, , [projection of A/ along the x-axis] can be obtained by taking cosine 
of the chord length, A l c [eq. (C-2)] as: 


A L, = A L cos 6 : _ 


( n 

oW • M f. 

0 Al ) 


= 2 R sin — = cos t 


l 2 J 

2 R \ 

' 2 R J 


Al c 



(C3) 


(i = l,2,3,....,«) 


Substitution of equations (Cl) and (C2) into equation (C3) yields the deformed domain length, A L t , 
along the deformed x-axis as equation (C4): 


A L =2 


+ C _i 

\ £ i-l +£ ij 


sin 


A / 

. 2 c i-l + 


COS 


0i-l + 


A / £ m + £. 
2 c t _ x + c ; 


(C4) 


(i = 1,2,3,.. ..,«) 


which is equation (24) in the main text. 
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APPENDIX D 


DERIVATIONS OF THE EULERIAN SLOPE EQUATIONS 

Equation (26) is reproduces as equation (Dl) in the following: 


r 

>y,-, 


dV 


[i+v 2 ] m 


< 


. x-x, , 
AL- 




*~*i- 

AL, 


dx ; (Xj- j < x < x f ) 


(Dl) 


On the left hand side of equation (Dl), the lower limit of the integration, V t _ x , is defined as the value of 
V at x M , namely, V t _ x = E(x M ) = tan 0(x.) = tan 0 t _ x . 

Integrations of both sides of equation (Dl) yield (refs.l, 17,): 

y K-i 


S i - 1 


(x-Xi-l) 


Vl+P 2 'fi+vli ( c i-i“ c i) 


AL, 


AL, 




(AL,) 2 

(*-Xi_i) 


• log, - X i -l) +C i-ij- l°ge C /-i| - 0 


+ AL, 


( c ii- c i) 

(W-l -C/fi-l) 


( C »-l ~ C i) 


log, 


(g/-l ~ C i) 
AL,. 


(x-X M ) + C,._ 


~/-l 


(D2) 


(•Vi < X < X,.) 


After simplification, equation (D2) becomes equation (D3): 


V = sf\ + V 


2 I P, 


i - 1 




C M “ C i 


-(x-x m ) + AL ; - 


(«h -^) 2 


-log, 


A A 

-3--1 

v c M y 


x -^ + i 


AL, 

(*;-i ^ x < x,.) 


(D3) 
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Squaring both sides and rearranging, and then taking the square root, the slope, V(x) = tan#(x) can be 
expressed as: 


M 


tan 6^ 


tan#(x) = 


ijl + tan 2 Q i 




-1 


C z-1 “ C z 


( x - x M ) + AL. 




( C i - 1 ~ C i) 


log, 


2 *=>e 


f \ 

*--l 

v c m y 


x - X, 


AL, 


izL + l 


(D4) 


tan 


7 = = + ~ 

yi + tan 2 0 t _ , c ; 


*/-i - 


(x-x M ) + AL, 


g ,--i c ,- ~ g<c M 
( c /-i - c i ) 2 


log £ 


C; 


v c m y 


x - X 


1-1 


AL, 


+ 1 


M 


(X;_ 1 < X < X, ) 


At the deformed strain-sensing station, x = x t , we have ( x, - x w ) = A/., , and equation (D4) becomes 
equation (D5): 


tan 6, 


i - 1 


tan <9 = 


•/l + tan 2 


+ AL, 


£;_.c, -Ac, , , C, 

— L + 7^ T^log, — 

C z ( C z-1 C i ) C Z-1 


(D5) 


1- 


tan 6, 


i-i 




: + AL, 


+ tan 2 # 


/-I 


£ i - 1 £ i + £i-l£l g i c i-l | Q g c i 


Cm-C/ (c M -c ; ) 


C, 


Z-l 


(/ = 1 , 2 , 3 ,. 


which is identical to the Eulerian slope equation (29) for large deformations of nonuniform beams (uniform 
beams excluded). 


Log-Expanded Slope Equation 

Note from equation (C4) that, the logarithmic terms and the denominators containing factor, (c i - c t _ x ) , 
will cause mathematical indeterminacy (that is, 0/0) for the uniform beam case (c t = c M = c ) . The way to 
bypass such a mathematical indeterminacy problem is to expand logarithmic function in series form in terms 
of (c- - c i _ ] ) . This approach enables the canceling out of the factor, (c- - c t _ x ) , in the denominator of 

x — X- 

equation (D4). Keep in mind that in equation (D4), the two conditions, 0 < — < 1 and 

AL, 


f \ 

-S--1 

< 0 , caused the required condition, - 1 < 

7 x 

-s-i 

x - X,-i 

V C z-l y 


_V C <-1 J 

AL,. 


- 1 < 


expansion of the logarithmic term, log e 


f \ 

\ c i - 1 y 


x-x w+ i 


AL ; 


< 1 , to hold for series 


expressed in terms of (c i - c M ) (ref 17) 


as: 
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log, 


f \ 

C , 


V c i - 1 ^ al z 


X - X; 


i-1 


+ 1 


A A 

^-1 

v c m y 


X - X._! 1 

AL ; 2 


f \ 

-s-i 

v c m y 


x - X 


\2 


z-1 


AL ; , 


C, - C ; _, x - X,_| 1 

/■ \ 
C/-CM 

2 

/ \ 
x-x w 

c i- 1 AZ, 2 

V C ;-1 y 


l w J 


+ .... 




In light of equation (D6), the term Min equation (D4) can be written as equation (D7): 


M v \ i \ j £ i-\ C i £ i C i - 1 


-(x-x M ) + AL ; - 


(cm -c,-r 


log, 


-S-- 1 

V c ;-i ) 


x - x , 


M 


AL, 


+ 1 


-(x — x M ) + AL ( ■ 


(Cm -c,) 2 


C; - C M X - X M 1 

/ 'N 

C/ 

2 

f X - 'l 

2 

Cm AL ; . 2 

V C M y 


l J 



Log-expanded terms 


= - x_, ) - ^ - **- (x - jc,_, ) — ^ ^ (J - J r ’ } + .... 


c m ( c m “ c i) 

(x X -_! ) C,- - £ 

( S i-l c i-l ~ £ i C i-\ ~ £ i-\ C i + £ i C i -\ ) - 2 


Cm (Cm - c ;) 

(•* - A-i ) 


2Ai, 

(x-x M ) 2 

2AL : 


+ .... 


Lm( C M -c,)]- j 


-i-1 


(X~X,._!) 2 

2AL 


c m( c m -c,) 

After canceling out the factor, (c-_ x - c.) , the term Mbecomes: 

(x-x M ) 2 


+ .... 


t. £ i — 1 / \ £ i—\ £ i £ i £ i — 1 

M = — — (x - x M ) — — LJ - L 

c m c M 


2AL ; 


+ .... 


In view of equation (D8), the slope equation (D4) takes on the following form: 

tan 

•^/l + tan 2 # ; _i 


tan 6(x) ■ 


tan 0 M e ,_ : i . «mC, -fiC,., (^-^m) 2 

c„ ( m) " 2Ai, 


1 


tan 0 t _ x , v , CMC,-fi,c M (x-x M ) 2 


1 + tan 2 _ 


+ — (x-Xm)' 


"M 


C, 


M 


2AZ,, 


(D6) 


(D7) 


(D8) 


(D9) 


(x M ^x<x t ) 
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At the deformed strain-sensing station, x = x i9 we have (x z - - jt M ) = A L t , and equation (D9) becomes: 


M 


tan ft M + g M _ g,_|C, - g,c M AL,. 


tan ft 


•^/l + tan 2 ft_. 


'i-l 


'i-l 


V 


tan ft_! ft-,ft -£,c,_ t AL, 


1 + tan 2 0. 


-l 


c- 


i-i 


(DIO) 


(/ = 1 , 2 , 3 , ....,«) 


The term M in equation (DIO) can be further simplified as: 

ft-ift-ftft-i A L, 


M=^-AL : -- 
ft- 1 


"i-l 


ft-i ai . AL > ft-ift~ftft- 1 


-i-i 


= AL 
AL, 


^ 2 ft-ift-i “ft-ift + ftft-f 


2c 


M 


Yf \ 


- C r 


2 l — 

ft- 1 + ft 

_V C i - 1 J 



2 ft-i 


In view of equation (Dll), equation (DIO) takes on the final form as: 


tan ft 


i-l 


tan ft = 


1- 


tanft 


al, r 

f 

9 _ 

.A.1 

2ft- 1 l_ 


ft-J ' 

+ M ' 

Ef, ft 

2ft-! 

V 

ft-i- 


(Dll) 


(D12) 


(/ = 1 , 2 , 3 ,. ...,«) 


Equation (D12) is the log-expanded Eulerian slope equation (34) for large deformations of nonuniform 
beams including uniform beams as special cases. 
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APPENDIX E 


DERIVATION OF THE L AGRAN GIAN SLOPE EQUATION 

Equation (42) is duplicated as equation (El) in the following: 


| 'V dV _ r* £|1 ^ ,1 Al 

■E. J-\-V 2 J'H , " 


jc-x 


; (i M < * < x,-) (El) 


A/ 


In equation (El), V = dyl dx = tan d(x ) , and the lower limit of integration, V j , , on the left hand side is the 
value of V at the strain-sensing station, x = x M [that is, V t _ x = (dy I dx) t _ x = tan #(x ,_, ) = tan 6 t _ x ]. 
Integrations of both sides of equation (El) yield (refs. 1, 17): 


sin 1 V - sin 1 V t _ x 


(x-x M ) 

( C M “ C i) 

c M (f M -f,0 £ M (c M -c ; ) 


+ - 


A/ 


A/ 


(cm-c ,-) 2 

(A/) 2 


log, 


( C <-1 ~ C ,) 

M 


(x-x,_i) + c,_. 


- log, C ,-l 


g ' ; 1 ^ - (x - x M ) + A/ C, '- l(g, ~- 1 / <} gM 2 (C;M C,) log t 


(C/-1 -c,) 




( C ,-l -C,~) 

A/ 


(x-X M ) + C M 


"z -1 


(E2) 


After simplifications, equation (E2) becomes: 


sin -1 V = sin -1 V t _ x + g, \ (x - x M ) + A/%1^— ^ip-log e 


( c w ~ c i ) 

Equation (E3) is equation (43) in the text. 


( C i-1 - C i) 


( C i-1 - C i) 


A/c, 


Z-l 


(^M <X<X { ) 


(x - X,_, ) + I 


0,-1 < X < X,.) 


(E3) 
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